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SUMMARY 


I A  Holzer-Myklestad  type  of  procedure,  using  a  matrix  formulation,  is  de¬ 
veloped  for  the  determination  of  the  natural  vibration  characteristics  of  a 
pretwisted  rotating  blade  in  coupled  bending  and  torsion.  The  nonrotating 
blade  is  considered  as  a  special  case.  Results  of  a  limited  parametric  study 
are  presented^]  It  is  found  that  in  the  case  of  the  rotating  blade  there  can 
be  an  appreciable  effect  of  centrifugal  forces  in  coupling  the  bending  and 
torsional  vibrations. 


'In  order  to  investigate  the  effects  of  Coriolis  forces  and  the  nonlinear 
effects  of  large  angular  displacements,  a  study  is  made  on  the  basis  of  a 
simple  model]  Numerical  results  indicate  that  the  Coriolis  forces  may  intro¬ 
duce  sub s t ant i^^pSfs e  differences  between  bending  and  torsional  vibration. 
Limited  numerical  results  on  the  nonlinear  effects  indicate  that  these  effects 
decrease  slightly  the  frequency  of  the  characteristic  motions  as  determined 
from  a  linearized  analysis  and  introduce  some  coupling  between  the  character¬ 
istic  motions. 
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1.  INTRODUCTION 


In  a  previous  report^  the  natural  vibration  characteristics  of  rotating 
twisted  blades  were  studied  for  the  special  case  of  coincident  mass  and  elas¬ 
tic  axes.  This  eliminates  coupling  between  bending  and  torsional  vibration, 
and  the  problem  was  studied  as  one  in  bending  vibration  only.  Bending  de¬ 
formation  about  both  principal  axes  of  the  cross  section  was  considered. 

The  present  work  represents  an  extension  of  this  previous  work  to  the 
case  of  noncoincident  mass  and  elastic  axes,  that  is,  the  case  of  coupled 
bending  and  torsion..  This  case  has  already  been  treated  analytically  in 
rather  complete  fashion  in  Ref.  2,  the  problem  being  formulated  in  terms  of 
governing  differential  equations  and  also  in  terms  of  energy  principles.  How¬ 
ever,  very  few  results  are  presented  in  that  reference,  and  they  are  for  a  few 
special  cases  of  a  rather  restrictive  nature. 

In  the  present  work  a  different  analytical  approach  has  been  used.  It 
involves  essentially  an  extension  of  the  Holzer-Myklestad  method  for  deter¬ 
mining  the  bending  vibrational  characteristics  of  a  beam  to  the  case  at  hand. 
The  Holzer-Myklestad  method  had  previously  been  extended  by  Targoff^  to  the 
case  of  bending  of  twisted  rotating  blades  and  applied  in  Ref.  1.  It  was 
found  to  be  particularly  well -suited  to  automatic  digital  computation,  and, 
for  that  reason,  has  been  extended  in  the  present  work  to  include  torsion  as 
well,  and  has  been  applied  in  a  limited  parametric  study. 

An  effect  of  centrifugal  forces  in  coupling  bending  and  torsional  vibra¬ 
tion,  considered  initially  in  Ref.  2,  is  taken  into  account  in  the  present 
work.  It  arises  when  the  mass  and  elastic  axes  of  the  blade  are  not  coinci¬ 
dent. 


It  should  be  remarked  that  the  inclusion  of  torsional  deformation  compli¬ 
cates  the  effects  of  pretwist  and  rotation  considerably.  There  may  be  a  siz¬ 
able  steady-state  or  "pseudo  static"  torsional  deformation  of  the  rotating 
blade  in  some  cases.  This  is  due  to  centrifugal  twisting  moment  which,  in  the 
case  of  negative  pretwist  and  positive  pitch,  tends  to  twist  the  blade  nega¬ 
tively,  and  also  to  the  twisting  moment  associated  with  tensile  stress  in  the 
longitudinal  fibers,  the  so-called  "centrifugal  untwisting  moment."  These 
two  effects  oppose  each  other  in  the  normal  case,  and  the  extent  to  which  one 
or  the  other  predominates  depends  primarily  upon  the  amount  of  pretwist  and 
the  pitch  setting  of  the  blade.  An  analysis  of  this  deformation  and  presen¬ 
tation  of  some  results  are  given  in  Ref.  4. 

Additional  effects  relate  to  a  departure  of  the  torsional  stiffness  from 
the  value  provided  by  Saint  Venant  theory.  This  departure  is  associated  with 
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inclination  of  the  longitudinal  filers  of  the  blade  with  respect  to  the  elas¬ 
tic  axis,  due  to  both  pretwist  and  torsional  deformation.  The  normal  stresses 
in  these  fibers  can  be  seen  to  have  components  in  the  plane  of  a  cross  section 
and  to  exert  a  torsional  moment  about  the  elastic  axis.  They  arise  from  two 
sources.  Firstly,  there  are  normal  stresses  associated  directly  with  torsional 
deformation  that  are  present  even  in  a  nonrotating  blade.  These  stresses  may 
introduce  a  substantial  nonlinearity  into  the  torsional  stiffness .5 “5  Secondly, 
there  are  normal  stresses  associated  with  centrifugal  forces,  contributing  to 
the  torsional  stiffness  in  a  manner  which  is  essentially  linear  for  practical 
deformations:  that  is,  there  is  a  linear  relationship  between  torque  and  elas¬ 
tic  twist. 3;^  Some  theoretical  results  for  the  case  of  torsional  vibration, 
with  some  or  all  of  these  effects  included,  are  presented  in  Refs.  4  and  6. 

Because  of  the  possibility  of  substantial  pseudo-static  torsional  deforma¬ 
tion  and  nonlinearity  in  the  torsional  stiffness,  an  accurate  determination  of 
the  natural  frequencies  of  vibration  of  a  twisted  blade  should  be  based  on 
linearization  with  respect  to  the  pseudo-static  deformation.  This  has  not  been 
done  explicitly  in  generating  the  results  presented  in  the  present  report.  The 
values  of  pretwist  selected  must  be  interpreted  to  include  pseudo-static  de¬ 
formation.  This  facilitates  comparison  with  the  results  of  Ref.  1,  where 
pseudo -static  torsional  deformation  would  have  an  influence  on  bending  vibra¬ 
tional  characteristics,  and  where  the  values  of  pretwist  must  be  similarly  in¬ 
terpreted  to  include  such  deformation. 

Another  interesting  aspect  of  the  rotating  blade  vibration  problem  is  dis¬ 
cussed  in  Ref.  7*  It  is  shown  that  Coriolis  forces,  or  so-called  ’’secondary 
inertia"  forces,  associated  with  the  combined  vibrational  and  rotational  mo¬ 
tion  introduce  a  phase  difference  between  the  bending  and  torsional  vibration. 

In  order  to  investigate  this  effect  more  fully  and  to  investigate  the  non¬ 
linear  effects  of  substantial  angular  displacements  on  the  dynamic  character¬ 
istics  of  a  rotating  blade,  an  additional  study,  reported  in  Section  was 
conducted  on  the  basis  of  a  simple  model.  The  nonlinear  effects  considered 
are  those  associated  with  inertia  forces.  Nonlinearity  in  the  torsional 
stiffness,  as  discussed  above,  and  the  effects  of  centrifugal  tension  on  the 
pseudo -static  deformation  and  on  torsional  stiffness  are  not  included,  al¬ 
though  they  could,  in  any  extension  of  the  present  work,  be  included  without 
undue  complication. 
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2.  BLADE  ANALYSIS 


SYMBOLS 

A  =  GJe  +  Tk|  +  EBi(P' )2 
A  =  A/EI1q 

Bi,  B2  section  constants  defined  in  Appendix  A 

C  =  EB2p!/A  =  C 


E  Young’s  modulus 

Eli,  EI2  "bending  stiffness  at out  major  and  minor  principal  centroidal  axes, 
respectively 

Eli  =  EIi/EI !q ,  Els  =  Ejl^/EIsq 

e  distance  Between  mass  and  elastic  axis,  positive  when  mass  axis 

lies  ahead 

e  =  e/R 

e^  distance  "between  area  centroid  of  tensile  member  and  elastic  axis, 

positive  when  centroid  lies  ahead 

eA  =  eA/R 

eQ  distance  at  root  "between  elastic  axis  and  axis  about  which  Blade  is 

rotating,  positive  when  elastic  axis  lies  ahead 

eo  “  e0/R 

GJe  effective  torsional  rigidity 

GJe  =  GJe/EIlQ 

I^,  I  mass  moment  of  inertia  of  cross  section  about  £  and  Y)  axes,  re¬ 
spectively,  defined  so  that  corresponding  moments  for  an  element  dx 
are  I^dx  and  I^dx 

H  =  yp0R2  y  =  iT]/p0r2 
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kA 

polar  radius  of  gyration  of  cross-sectional  area  effective  in 
carrying  tensile  stresses  about  elastic  axis 

II 

I** 

kp, 

mass  radii  of  gyration  about  £  and  t]  axes,  respectively 

l 

length  of  blade  segment 

1 

1/ R 

Mi,  M2 

bending  moment  about  major  and  minor  principal  axes  of  cross 
section,  respectively,  when  centrifugal  tension  is  assumed  to 
act  along  undeformed  position  of  elastic  axis 

m 

mass  of  blade  segment 

• 

px,  p5,  Pn 

resultant  loadings  per  unit  length  in  the  x,£,q  directions, 
respectively 

Q 

resultant  torque  about  elastic  axis  at  any  cross  section 

^x’  %’  'It] 

resultant  torsional  loadings  per  unit  length  about  the  x,£,r] 
axes,  respectively 

R 

blade  radius 

T 

2 

centrifugal  tension,  ft  Ti 

N 

t= I p 

i^ixi 

i=l 

N 

Ti  - 1 p 

i*ixi 

i=l 

u 

displacement  in  the  x  direction 

Vi,  v2 

shearing  forces  in  the  direction  of  the  minor  and  major  principal 
axes  of  the  cross  section,  respectively 

x,y,z 

coordinate  system  which  rotates  with  blade  (Fig.  2.2) 

X 

II 

\x 

Y  =  (EB23’ 

1  )2/E Is 

b 

Y  =  Y/EIl0 


angle  between  major  principal  axis  of  cross  section  and  plane  of 
rotation^  either  in  the  undeformed  or  pseudo-static  state 


p'  =*  dp/dx 


P1  =  p'R 

Ap  increment  in  p  "between  hlade  segments 

72  =  ei1o/ei2o 

5  9  5  displacements  of  the  elastic  axis  in  the  y  and  z  directions; 

^  respectively 

&x;  52  displacements  of  the  elastic  axis  in  the  direction  of  the  minor 

and  major  principal  axes  of  the  cross  section;  respectively 


T) 


© 


coordinates  in  direction  of  minor  and  major  principal  axes, 
respectively 

total  twist  in  hlade  between  x  =  0  and  x  =  R,  ©  =  -Rp1 


\  =  0)n/p0R4/EIi0 

=  a  n/  p0r4/eIio 

p  mass  per  unit  length  of  hlade 


p  =  p/p0 


* 

O) 

n 

t  ] 
{  ] 


torsional  displacement;  positive  when  leading  edge  is  up 

natural  frequency  of  hlade  vibration 

rotational  velocity 

rectangular  matrix 

column  matrix 

Other  symbols  are  defined  in  the  text . 


Subscripts 

n  order  of  natural  mode 
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0 


value  at  x  =  0 


T  value  at  x  =  R 

(  )’,  (  )M  differentiation  with  respect  to  x,  except  in  Appendix  C. 

BASIC  MATRICES 

The  governing  differential  equations  of  motion  for  a  rotating  blade  with 
offset  mass  and  elastic  axes  have  been  derived  and  are  reported  in  Ref.  2. 

These  equations  are  repeated  in  Appendix  A.  In  the  present  report  these  equa¬ 
tions  have  been  adapted  to  a  matrix  formulation  which  permits  rapid  numerical 
analysis.  This  method  is  essentially  an  extension  of  the  one  presented  in  Ref. 
1. 


The  coordinate  axes  of  the  blade  are  shown  in  Fig.  2.1.  The  cross  section 
coordinates  and  displacements  are  shown  in  Fig.  2.2.  The  blade  is  divided  into 
a  number  of  spanwise  segments,  not  necessarily  equal  in  length.  The  mass  of 
each  segment  is  assumed  concentrated  at  its  center,  and  the  bending  stiffnesses, 
Eli  and  EI2,  the  torsional  stiffness,  GJe,  and  the  angle  of  incidence,  p,  are 
assumed  constant  between  masses,  appropriate  average  values  being  selected. 

The  built-in  twist  is  accounted  for  by  relative  rotations  of  adjacent  uniform 
bays  (between  masses)  about  a  spanwise  axis,  the  change  in  angle  Ap  being 
equal  to  the  total  twist  in  a  segment  and  occurring  just  outboard  of  the  mass 
(Fig.  2.3). 

The  quantities  Vi,  Mi,  Si,  Si,  V2,  M2,  62,  82,  Q,  and  <f>  (Fig.  2.2),  which 
apply  when  the  beam  is  at  its  maximum  displacement  in  a  free  vibration,  are 
defined  at  stations  along  the  beam  and  may  be  represented  at  any  station  in 
the  form  of  a  column  matrix: 


CM 


=  < 


rv£l 

Mi 

Bi 

Si 

V2 

M2 

82 

82 

Q 


(2.1) 


The  elements  of  this  matrix  will  vary  along  the  beam  in  such  a  manner  that 
the  variation  can  be  considered  to  occur  in  a  series  of  steps.  Moving  from  the 
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Fig.  2.2.  Nomenclature  and  sign  convention  for  cross-section  coordinates 
displacements ,  tending  moments,  and  shears. 


Fig.  2.3-  Blade  segment  rotation. 


tip  toward  the  root  of  the  beam,  the  change  in  [A]  occurring  from  a  station 
immediately  outboard  of  one  mass  to  a  station  immediately  outboard  of  the  next 
mass  can  be  broken  down  into  three  steps,  the  first  involving  movement  across 
the  mass,  the  second  involving  movement  from  one  end  to  the  other  of  a  weight¬ 
less  uniform  bay,  and  the  third  involving  movement  across  the  discontinuity  in 

P. 


The  relationship  between  the  [A}  matrices  as  they  apply  at  the  two  ex¬ 
tremes  of  this  travel  can  be  represented  as  follows: 


(A}n+1  [R][E][F][A}n 


(2.2) 


where  [F],  [E],  and  [R  ]  are  rectangular  matrices  representing  linear  rela¬ 
tionships  corresponding  to  the  three  steps  discussed  previously. 


The  [F]  matrix,  relating  the  (A)  matrices  on  either  side  of  a  concen¬ 
trated  mass,  is  written  as  follows: 


F 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0  0  F  i4 

1  F23  F24 

0  10 
0  0  1 

0  0  f54 

0  0  0 

0  0  0 

0  0  0 

0  0  F94 

0  0  0 


0  0  0 

0  0  0 

0  0  0 

0  0  0 

10  0 
0  1  F67 

0  0  1 

0  0  0 

0  0  0 

0  0  0 


F  i8  0  F110 

0  0  0 

0  0  0 

0  0  0 

F58  0  F510 

F68  0  0 

0  0  0 

10  0 
F98  1  F910 

0  0  1 


(2.3) 


where 


F14 
F 18 
F 1 10 
F23 


ip(a>2+ft2sin2p) 

2 

-  ipft  sin  p  cos  p 

ipe  {oo2+n2(sin2P-cos2p) )  -  ipeQ^cos  p 


(2.4) 
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F24 

= 

-  pixft2 

F210 

= 

2 

-piexft 

r10 

= 

F 18 

Fss 

= 

pi(tD2+fi2C0S2p) 

F510 

= 

-  2p^eH2sin  (3  cos  p  -  pi^eQ  sin  (3 

^67 

= 

^68 

= 

F24 

F94 

= 

2  2 

pie(o)  sin  p) 

F98 

= 

-  pief^sin  (3  cos  (3 

O 

H 

0> 

P»H 

= 

)  (cos2 (3-sin2 (3)  ffi2-ipee0fi2cos  (3 

The  derivation  of  the  elements  of  this  matrix  is  given  in  detail  in  Ap¬ 
pendix  B,  except  for  the  contribution  of  centrifugal  force  coupling,  which  is 
treated  separately  in  Appendix  D.  It  is  seen  that  only  the  shear  forces,  bend¬ 
ing  moments,  and  torque  are  changed,  since  there  are  no  discontinuities  in 
slope  or  displacement.  The  changes  in  shear  force  are  due  partly  to  the  in¬ 
ertia  force  associated  with  the  vibrational  motion  of  the  mass  and  partly  to 
the  component  of  centrifugal  force  normal  to  the  undeformed  position  of  the 
elastic  axis.  Part  of  the  change  in  torque  is  related  to  the  change  in  shear, 
since  the  mass  and  elastic  axes  do  not  coincide,  and  part  is  due  to  the  inertia 
force  associated  with  the  torsional  vibrational  motion.  The  change  in  bending 
moment,  except  that  associated  with  centrifugal  force  coupling,  is  fictitious 
and  arises  from  a  special  feature  of  the  analysis.  This  feature  involves  the 
replacement  of  the  component  of  the  centrifugal  force  parallel  to  the  unde¬ 
formed  position  of  the  elastic  axis  by  an  equal  force  along  the  line  of  the 

undeformed  axis  and  an  appropriate  couple  to  provide  static  equivalence.  The 

changes  in  bending  moment  indicated  in  the  [F ]  matrix  are  then  due  only  to  the 
applied  couple,  the  moment  due  to  the  force  applied  along  the  undeformed  axis 

being  accounted  for  in  the  [E ]  matrix.  When  moments  due  to  both  sources  are 

considered,  the  discontinuity  in  bending  moment  disappears.  Note  that,  on 
the  basis  of  this  procedure,  the  bending  moment  at  any  station  is  not  M,  but 
rather  M  plus  the  moment  of  the  tensile  force  T  acting  along  the  undeformed 
elastic  axis. 

The  elements  in  the  [E  ]  matrix  are  found  by  the  solution  of  the  differ¬ 
ential  equations  of  combined  bending  and  torsion  of  the  weightless  uniform 
bay  between  masses.  These  equations  and  their  solutions  are  given  in  Appen- 
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dix  C.  The  resulting  [E  ]  matrix  is: 


1 

0 

0 

0 

0 

e21 

1 

0 

0 

0 

e31 

E32 

E33 

E34 

E35 

E41 

E42 

E4  3 

E44 

E45 

0 

0 

0 

0 

1 

0 

0 

0 

0 

E65 

e71 

E72 

e73 

E74 

E75 

E8i 

E82 

E83 

E84 

E85 

0 

0 

0 

0 

0 

Eioi 

e102 

E103 

El04 

E105 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

E36 

E37 

E38 

E39 

E310 

E46 

E47 

E4  8 

E4  9 

E410 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

E76 

E77 

E78 

CD 

0 

E86 

E87 

E88 

E89 

0 

0 

0 

0 

1 

0 

E106 

E107 

E108 

E109 

E1010 

where ,  if  we  define 

Pl  = 

(ai-a3) (pf-pf ) 

P2  =  (p!-ai)(pl-ag) 
(ai-a3)(pi-p2) 
2  2 

p  =  (Pi-aa)(pg-ai) 
(ai-a3 )  (pf-p! ) 

p  _  (p|-a3)(pl-a3) 

(a1-a3)(p?-p|) 


(2.6) 


the  components  of  E  are  given  below.  The  quantities  a-j_,  and  fj_  are  de¬ 
fined  in  Appendix  C. 


E2 


E31  =  -  T^r- (a2a6P4 -a3aaP3 )cosh  Pll  +  [-a2a6P4+a3a8(P3+l)  ]cosh  p2i-a3aej- 


fgEIi T 


E32  =  +•  — sinh  Pii  +  — sinh  p2i 

EI1P1  EI1P2 


(2.7) 


E33  =  -  P3  cosh  pi i  +  P2  cosh  P2 l 
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E34  =—  aiEIiE32 


E35  =  - 


a£P4 


a3a4EI2  \pi 


1  1 

sinh  pii  -  : —  sinh  p2i 


P2 


E36  =  (cosh  Pl^-cosh  p2l) 


E37 


a3EI2E35 


E38  =  “  a3EI2E38 


E39  - - ,a|  g,  [~-(pi-a3)pf  cosh  pii+(p|-a3)pf  cosh  p2i]  +  ^3- 

f2(pf-P 1)  L  J  fa 

E310  =  -  asEIxEas  +  E36W^^ 


E4i 


—  {(-c 

f2El!  L 


-a3a8P3+a2a6P4 )  —  sinh  pii 
Pi 

ta3a8(P3+l)-asa6P4 ]  ““  sinh  p2i-a3a8i 


E42  = 


1 

Eli 


22  2 

-  Pi(p2-a3)+a3(p2-a1) 


(a.i-a3)pip| 


Pi 


cosh  Pii  +  ^§-  cosh  p2i 
P2 


E43  =  -  EIiE^ 


2  2 


E44  =  -  a-3 ■—  -  ai^,a  cosh  Pii  +  ai^g'  cosh  p2i 


22 

P1P2 


Pi 


^2 

P2 


E45  =  - a2|42„T"  ■  Rpi-pi)  +  picosh  Pii  -  p-fcosh  p2i] 

a3a4pip2EI2  L  J 


E4b  -  E35 


E47  =  -a3EI2E45 


E48 


E37 


(2.7) 
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^49 


a2 

f2(pi-pi) 


2  P  2 
£2  sinh  Pii+(p2-a3)  Ei  sinh  p 2i 
Pi  P2 


+  a2a3l 

f2 


E410  =  a2EIiE42  +  E46w(N) 


Ess  -  ^ 


E71  =  *  a'  YEI~  [(a2a6p2-a3a8Pi  )Pi  sinh  pii  +  (-a2a6P3+a3a8Pi)p2  sinh  p2i 


E72  =  +  a4Pl  (cosh  pii-cosh  P2!) 

a2EIi 


E73  =  (-Pi  sinh  Pii+P2  sinh  p2i) 

a2 


E74  =-aiEIiE72 


E75 


E76 


e77 


='^  (E77_1) 

- - i=-  (P2P1  sinh  Pii-P3P2  sinh  p2i) 

a3EI2 

=  P2  cosh  pii  -  P3  cosh  p2i 


(2.7) 


E78  =  "  a3El2E73 

E79  =-  r -(pf-ai)  sinh  pi*  +  (pf-ai)  sinh  p2il 

fa  ( Pi -Pa)  L  Pa  Pi  -I 

(N) 

E7io  =  -^a2EIiE72  +  E76W 

E31  =  — — -  I (a2aeP2-a3a8pi )cosh  pif  +  ( -a2a6P3+a3a8pi )cosh  P2^“a2a6 

a2f2EIi  L 

Es2  =  -  g&T  (fe sinh  Pli  ■  h sinh  P2*) 

Eq3  =  ~EIiE72 
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E84  =  aiEIiE82 


‘85 


a3EI2 


(Ee7-i) 


E8s  -  a3EI2  (Ess-l) 


E8t  =  —  —  sinh  pii  +  sinh  p 2i 
Pi  Pa 


E88  =  P2  cosh  Pi£  -  P3  cosh  p2i 


E89  =  — — 27  r-(p!-ai)pl  cosh  Pii+(pl-ai)Pi  cosh  pail  + 

fa(Pi-P2)  L  J  2 


E8io  =  a2El!E82  +  E86w 


(N) 


E 


101  = 


a2f2EI 


_  |~  (pi (a2a8P2-a3a8Pi )  sinh  pii 

i  L  pi 


®102 


^P^~.  ( a2a8P3 ~a3a8Pi )  sinh  p2^“^a2a3ag^ 

Do  J 

-  [a3(pi-pl)-(pi-a3)p!cosh  p1i+(pl-a3)pfcosh  p2*J 


P2 

Pi 


a2pip2EI; 


r~  2  P  ** 

Eio3  =  -  —  [-  ^Pl~as-  sinh  pif  +  ^Pa~a3^  sinh  p2i 

a2  L  Pi  P2  J 


E104  =  aiEIiEi02 


ElOS  " 


Eios  = 


a3a4pip2EI2 


a3a4EI2 


|j-ai(pi“p2  )  +  (pi-ai  )p2  cosh  pii-(p2~ai  )pi  cosh  p2 
\ sinh  Pii  -  (pg"ai.I  sinh  p2£ 

L  Pi  P2 


Ei07  =  -a3EI2Eios 


Eio8  =  a3EI2EiQ6 


15 


i„»  -  k  ‘  -  ^4=^ 

(P1-P2JP1 


sinh  pii  + 


(pj-aj.  )(pj-a3)pi 

(p?-pi )P2 


E1010  ~  pp 


ai 


P1P2  (ai-a3) 


Pi ( pi “a3 ) “a3 ( pi ~ai 


+  pi  ! “(pl“a3 )P2  cosh  Pli 


sinh  x>sl\ 

(2-7 ; 


ll 


+  (P2-a3)pi  cosh  p2ii  f  +  Eioew 


ir 


,(N) 


I  t 

Note  that  Si  and  62  are  positive  for  increasing  deflection  in  the  positive  x 
direction. 


The  [R  ]  matrix  serves  to  rotate  the  coordinate  axes  through  the  angle  A0 
and  is  written  as  follows: 


[H] 


cos  A3 
0 
0 
0 

sin  A3 
0 
0 
0 
0 
0 


0 

cos  A3 
0 
0 
0 

sin  A3 
0 
0 
0 
0 


0 

0 

cos  A3 
0 
0 
0 

sin  A3 
0 
0 
0 


0 

0 

0 

cos  A3 
0 
0 
0 

sin  A3 
0 
0 


-sin  A3 
0 
0 
0 

cos  A3 
0 
0 
0 
0 
0 


0  0 

-sin  A3  0 

0  -sin  A3 
0  0 

0  0 

cos  A3  0 


0  cos  A3 

0  0 

0  0 

0  0 


0  0 

0  0 

0  0 

-sin  A3  0 

0  0 

0  0 

0  0 

cos  A3  0 

0  1 

0  0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


(2.8) 


METHOD  OF  SOLUTION 


By  a  successive  multiplication  of  the  appropriate  matrices,  a  linear  re¬ 
lationship  can  he  established  between  the  (A]  matrices  at  the  root  and  tip  of 
the  beam 


f^root  [C](A}^p 


(2.9) 


Recognizing  that  the  shears,  bending  moments,  and  torque  are  zero  at  the 
tip  of  the  beam,  the  (A}^j_p  matrix  can  be  reduced  to  a  five -element  matrix, 
and  the  corresponding  five  columns  can  be  eliminated  from  the  first  [F  ]  matrix 
at  the  tip  of  the  beam;  successive  multiplications  will  then  yield  a  10  x  5 
matrix  product. 

In  order  to  satisfy  the  boundary  conditions  at  the  root  of  the  beam,  the 
determinant  of  a  5  x  5  matrix  formed  from  appropriate  elements  of  the  [C  ] 
matrix  must  equal  zero.  For  example,  for  a  cantilever  blade  the  third,  fourth, 
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\ 

seventh,  eighth,  and  tenth  rows  form  the  5*5  determinant,  and  for  a  fully 
articulated  hlade  with  torsional  restraint  the  second,  fourth,  sixth,  eighth, 
and  tenth  rbws  form  the  determinant.  Other  boundary  conditions,  such  as 
elastic  restraint  at  the  root,  can  he  handled  easily. 

The  elements  of  this  determinant  will  he  polynomials  in  (x£ ,  and  upon 
expansion  a  polynomial  equation  in  oo2  will  he  obtained.  In  principle,  the 
natural  .frequencies  of  the  hlade  could  he  determined  by  solving  for  the  roots 
of  this  equation;  however  ,  such  a  procedure  is  far  too  cumbersome  to  he  feas¬ 
ible  . 

A  more  practical  procedure  involves  the  introduction  of  trial  values  of  o 
into  the  various  [F  ]  matrices  and  evaluating  the  elements  of  all  matrices 
numerically.  The  matrix  multiplications  can  then  he  carried  out  numerically, 
and  the  appropriate  determinant  evaluated.  The  value  of  this  determinant, 
which  may  he  termed  the  "residual, 11  may  then  he  plotted  versus  cd  or  <x£  and 
the  location  of  the  zeros  of  the  residual  will  determine  the  natural  fre¬ 
quencies  of  the  hlade. 


STEADY-STATE  DEFORMATION 

As  pointed  out  in  the  introduction,  there  may  he  a  sizable  steady-state 
or  "pseudo-static"  torsional  deformation  of  the  rotating  blade  in  some  cases. 
The  loadings  which  produce  this  deformation  are  given  in  Appendix  B  along 
with  those  induced  by  the  lateral  and  torsional  vibratory  motion.  It  is 
possible  to  determine  this  pseudo-static  deformation  and  to  then  find  the 
natural  frequencies  based  on  linearization  with  respect  to  the  pseudo-static 
deformation.  In  the  numerical  results  which  follow  this  has  not  been  done 
explicitly.  The  values  of  pretwist  selected  should  be  interpreted  to  include 
the  pseudo -static  deformation. 

In  order  to  determine  the  pseudo -static  deformation  let  us  define  the 
following  matrices: 


=  column  matrix  of  blade  variables  just  outboard  of  mass  i 

{^}jL  =  column  matrix  of  blade  variables  just  inboard  of  mass  i 

[F0].  =  matrix  [F]  with  cd  =  0,  across  mass  i 

[d}#  =  column  matrix  of  steady  state  quantities  across  mass  i 

(g).  =  column  matrix  of  steady  state  quantities  across  bay  between 

masses  i  and  i+1 

[D]^  =  [R  ] ^  [E  ] i  across  bay  between  masses  i  and  i+1. 
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Then  it  follows  that 


and 


Starting  at  the  root  where 


^^root  ^n+1  ’ 


we  have 


where 


[H]  ;  (h}1  UA).! 

lo]  ;  "i"  |  1 ' '  f  > 

—  L  vj 


'[H] 

(h) 

n 

-  If 

‘[DJi 

(g)i 

"[F0]  1  fd)± 

pj 

i=l 

;i?r 

i  -f" 

UQJ  1  ~i'~_ 

and  in  which  i  decreases  as  one  proceeds  from  left  to  right. 
Equation  (2.12)  may  he  written 


(2.12) 


tMroot  =  PHMi  +  (h) 


(2.12a) 


Satisfying  the  boundary  conditions  at  the  root  and  the  tip  of  the  blade,  Eq. 
(2.12)  may  be  reduced  to 


[hOO]  (a(D)  = 


-  (h(K)}  , 


(2.13) 


where 


tad)} 


and  where  K  =  1  corresponds  to  a  fixed  root,  and  K  =  2  corresponds  to  a  fully 
articulated  root  (Mx  =  Ms  =  0) .  [H^K) ]  is  a  square  matrix  of  order  5,  and 

(h(K)}  Is  a  five  element  column  matrix. 

In  the  case  of  a  fixed  root,  [h(D]  is  obtained  hy  deleting  rows  1,  2,  5> 
6,  and  9  and  columns  1,  2,  5,  6,  and  9  from  [H],  and  (h^1))  is  obtained  by 
deleting  rows  1,  2,  6,  and  9  from  (h).  Similarly,  in  the  case  of  a  fully 

articulated  blade,  [h(2)_]  is  obtained  by  deleting  rows  1,  5>  1>  9  and 

columns  1,  2,  5,  6,  and  9  from  [H ] ,  and  (h(2)]is  obtained  by  deleting  rows  1, 
3*  5;  1)  and  9  from  [h}. 

Equation  (2.13)  may  be  solved  for  (A^1)},  and  (A]  then  determined  for  all 
stations  by  applying  Eqs.(2.10)  and  (2.11),  starting  at  the  tip  and  progress¬ 
ing  toward  the  root. 

The  matrices  (d)  and  (g]  are  each  ten-element  column  matrices  which  are 
obtained  from  the  steady  state  terms  in  Appendices  B  and  C.  From  Appendix  B, 


r 

61 

Bi 

e>2 

62 

( 


(2.14) 


di 

=  -  P  1$ 

2  sin  g(e0- 

+e  cos  p) 

a 2 

II 

cn 

11 

=  d7  =  d8 

=  dio  =  0 

d-5 

=  pift2  ' 

cos  p(eQ+e 

cos  p) 

de 

=  -pi  X 

etf 

d8 

=  (V1 

^)i-Q2  sin 

P  cos  p  -  pJLeeQQ?  sin  P  . 

Appendix  C  shows  that  the  (g)  matrix  can  be  derived  from  the  [E]  matrix  if 
the  terms  involving  Mg  and  Q  are  extracted  and  M2  and  Q  are  replaced  by  Te^ 
andTk^p*,  respectively.  Thus, 


Si 


Ei2TeA  +  Ei9T4P' 


(2.16) 
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NONDIMENSIOML  FORM 


It  is  convenient  and  desirable  to  treat  the  problem  in  nondimens ional 
form.  The  (A)  matrix  can  be  redefined  in  terms  of  nondimens ional  forces, 
moments,  and  deformations  as  follows: 


(A)  = 


ViR2/EI1o 
MiR/EI1q 

Bi/R 

v2r2/ei1o 

^M2R/EIi0 

t 

62/R 

QR/EIln 


4> 


> 


(2.17) 


The  corresponding  nondimens ional  form  for  the  [F]  matrix  follows: 


[F] 


0 

F23 

1 

1 

0 

0 

0 

0 

0 

0 


F14 

F24 

0 

0 

F54 

0 

0 

0 

F94 

0 


0 

0 

0 

0 

0 

^67 

1 

0 

0 

0 


Fie 

0 

0 

0 

F58 

Fee 

0 

1 

F98 

0 


Fno 

0 

0 

0 

F510 

0 

0 

0 

F910 

1 


(2.18) 


where 
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_  _ _ _  g  g  ^ 

f14  =  -£pU  +p  sin  0) 

F is  =  -  ipp2  sin  p  cos  p 

Fno  =  pie{X2+|i.2(sin2p-cos2p) )  -  pfe0n  cos  P 

f23  =  eg^+  ii^V8) 

f24  =  -  P*xp2 

F2io  =  -  p2piex 

F54  =  F  is 

(2.19) 

F58  =  p7(X2+M-2cos2P) 

F510  =  "  pip2sin  p(2e  cos  P+e"0) 

F67  =  +  7^(\2+n2) 

12  s 

^68  =  ^24 

_  —  _  2  2  2 

F94  =  pie(X  +1-1  sin  p) 

F9S  =  -  ple|i2sin  p  cos  p 

Fgio  =  (l^+lTj)iX2  +  M2!  cos  2&(ir]-It>)  -  pie0Ip2cos  p 


19 


And  the  corresponding  [E  ]  matrix 


1 

0 

0 

0 

0 

l—l 

01 

IH 

1 

0 

0 

0 

E3i 

E32 

E33 

e34 

E35 

E4i 

E42 

E4  3 

e44 

E45 

0 

0 

0 

0 

1 

0 

0 

0 

0 

E65 

E71 

E72 

E73 

e74 

E75 

E81 

E82 

CO 

CO 

m 

E84 

E85 

0 

0 

0 

0 

0 

E101 

E102 

E103 

E104 

Eiq5 

where 
E21  =  i 

E31  = 

E32  = 

E33  =  -  P3  cosh  pi  /  +  P2  coshTs^ 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

E36 

E37 

CD 

lw 

E39 

E310 

E46 

E47 

E48 

E49 

E410 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

E76 

E77 

e78 

E79 

E710 

E86 

E87 

Eb8 

Es9 

^810 

0 

0 

0 

1 

0 

E106 

Eio7 

E108 

El09 

1 

(2.20) 


4(a2a6P4-a3a8P3)cosh  p  f+  [-a2a6P4+a3a8(P3+l)  ]cosh  p2i-a3a8 
12EI 1  l 


sinh  pii  -  -?2_~  sinh  p2i 
EI1P1  EI1P2 


E34  =  aiEIiE32 


^35  ~ 

E36  = 

E37  = 

E38  = 

E39  = 


-  7  (zr~  sinh  pi l  -  i—  sinh  p2i  ) 

a3a4Sl2  \Pi  P2  J 


(2.21) 


Z-.,.g2^4 ■-  (cosh  Pli-cosh  p2i 
a3a4EI2 


-  ^2  e35 

♦V2 


^2  e36 

7 


w2 


- ag  -  f(p?-a3)pl 

T2(p5-p1)  L 


i|  cosh  pii-(pi-a3)pf  cosh  p2i  -  ^ 

J  fs 
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E3io  =  aaEIxEss  +  h2wWe 


36 


E4I 


(-a3a8P3+a2a6P4)  sinh  pi l 


- 

faEIi  t 

+  a3ai  ( P3+l )  -^aeP^-i^s  inh  p27-'5:3a8 1  r 


E42  =  — 

EIo.  L 


_2  ,_2  _  ,  _  ,_2  _ 


pl(pg-a3)+a3y-a!)  _  P3  cosh  j  +  P|  cogh  -  f 
(ai-a3)p?pl  pf  Pi  J 


E43  =  EI1E32 


a3(P!-p|,)Pi  .  SiPa  cosh  -1J  +  ii|2  cosh  52j 

Ijpl- 


^44 


—  2—2 
PlP2 


p? 


pi 


E. 


45  = 


a2P47 


a3a4pfpiEI2 


2  2  _  ^3  ^ 

p2)+P2  cosh  pii~pi  cosh  P2/ 


E4Q  =  -  E35 


E. 


47  ~ 


^E45 


E4s  = 

E49  “ 

E4I0  “ 


E37 


&2 

-  -2 

f2lPl-P2 


—  _ 2 

-(5f-a3)  Is-  sinh  Pii+CpI-'Ss )  if"  sinh  p2I 

L  pi  P2 


ll  a-ir,h  -„rl  +  a#37 

J  fa 


_ _  a_(N)_ 

IL2EI1E42  +  p  W  E40 


E, 


65 


E71  =  — 14_  r(-a2a6P2+a3a8Pi)pi  sinh  pii+(a2a6P3~a3a8Pi)p2  sinh  p2 

anf^KT,  L 


a2f2EXi 

E72  =  (cosh  'pii-cosh  p27) 

a3EIi 


E73  =  (-"Pi  sinh  ^ii+p2  sinh  p2*  ) 

3-2 


E74  =  aiEl!E72 


(2.21) 


<1 
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E?5 

e76 

E77 

e78 

E79 

E710 

Esi 

Es2 

E83 

Eq4 

E85 

^86 

E87 

E88 

E89  = 

E8iq  = 


~Z —  (I-E77  ) 

a8El2 


(-P2P1  sinh  Pii+p3P2  sinh  p2£) 

a3EI2 


P2  cosh  pii  -  P3  cosh  p2i 


as|l2  e76 
7 

Rpi-aiPsinh  Pii-(p2-a1f  1sinh  p2il 

f2(pl-pa)  L  -I 

-  -  p — (N ) — 

a2EIiE72  +  \i  W  E75 


,a4 


a2f2EIi 


L(a2a6P£ra3a8Pi)cosh  pii+(  -a2a6P3+a3a8Pi  )cosh  p2i-a2a6 


-*jLX' -  (=-  sinh  +  i-  sinh  p2i) 

a2EIi  \Ji  p2  J 


EIiE72 


{2.21 


a-)  Eli  E82 


-  _7i 

a3EI2 


(E87+i ) 


a3EI2 


(E88-l ) 


-  ~~  sinh  pii  +  ^2.  sinh  p2i 
Pi  P2 

P2  cosh  pii  -  P3  cosh  "p"2T 


- . ,ztr  ~ZZ~  [-(Pi-ai)pl  cosh  Pi/+(pf-ai)pf  cosh  p2i|  + 

f2(pi-p2)  L  J  f2 

_ -  2-(K)_ 

a2EIiE82  +  p  W  E86 
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rv-2 


E101  “ 


_2 


(pi-a3),  (a2aeP2-a3aspi)sinh  p if 


Pi 


(P2^a3 ),  (a2a6  P3-a3a8Pi )sinh  P2i 

P2 


] 


El 02  “ 


E103  - 


Pi 


aaPiplEI 


—  [a3(p?-pi)-(pf-'a3)pl  cosh  pii+(pl-a3)pf  cosh  p2f  +a2a3a6ij 
Ii  L 

rii  1  -  ^-3_1  airih  0-1  il 


Pi  (Pi”a3.)„  sinh  Pl£  -  .(Pg^-sJ.  sinji  Pli 

a2  L  Pi  P2  J 


E104  =  aiEIiEi02 

E105 


•106 


‘1q7 


=  - -  ■  ("-ai(pi“P2 )  +  (pi“ai )p2  cosh  Pi-Mpl-ai  )pi  cosh  p2lj 

a3a.4plp2EI2  (2.21) 

L_  Rii 
EI2  L  p, 


72pt. .  sinh  piF  -  4^  sin*  p2i 

a3a4EI2  L  Pi  P2  J 


aaE^Exos 
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The  nondimension  quantities  used  above  are  defined  as  follows: 
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where  N  denotes  the  number  of  masses  outboard  of  the  bay  under  consideration. 


The  (d)  matrix  has  the  components 


di  =  ip p2  sin  p(e0+e  cos  p) 
d2  =  d3  =  d4  =  d7  =  d8  =  di0  =  0 

d5  =  Ipp2  cos  p(e0+e  cos  p) 
d6  =  -  ipxep2 

dg  =  sin  p  cos  p  -  ipee0p2  sin  p 

and  the  components  of  [g]  are 

gj_  =  (^(Eisf^+Eisik^p*) 


(2.23) 


(2.24) 


For  the  case  of  zero  rotational  velocity,  the  [F]  matrix  is  obtained 
directly  by  substitution  of  \i  =  0.  When  this  substitution  is  made  in  the  [E] 
matrix,  some  of  the  elements  are  found  to  be  of  indeterminate  form  and  a 
limiting  process  must  be  applied.  This  results  in: 
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NUMERICAL  RESULTS 

A  program  of  computations  was  performed  for  two  representative  cantilever 
Llades.  The  properties  of  these  blades  are  given  in  Table  1.  They 
were  chosen  to  have  the  same  bending  properties  as  blades  for  which  numerical 
results  are  reported  in  Ref.  1.  For  these  two,  the  section  constants  Ii,  I2, 
Bi,  B2,  GJe,  k^  correspond  to  those  for  a  thin-walled  rectangular  section,  and 
it  was  assumed  that  some  nonstructural  mass  was  distributed  in  such  a  way  as 
to  provide  an  offset  between  the  mass  and  elastic  axes,  and  to  provide  suffi¬ 
cient  mass  moment  of  inertia  to  make  the  uncoupled  first  torsional  frequency 
and  the  second  uncoupled  flapwise  frequency  coincide.  The  result  in  both  cases 
is  a  lightly  coupled  system  as  far  as  flapwise  bending -torsion  is  concerned. 

In  both  cases  the  blades  were  divided  into  ten  segments,  the  cantilever 
root  condition  was  applied,  and  the  four  lowest  frequencies  were  determined. 

A  range  of  values  of  pretwist  and  rotational  velocity  was  chosen,  and  the  re¬ 
sults  are  presented  in  Figs.  2 .4-2 .6.  In  addition,  results  for  beam  No.  2 
with  rotary  inertia  neglected  are  presented  in  Fig.  2=6,  and  with  centrifugal 
force  coupling  neglected  in  Fig.  2.7° 


DISCUSSION  OF  RESULTS 

The  influence  of  twist  on  the  natural  frequencies  of  nonrotating  blades 
is  shown  in  Fig.  2.4.  The  fundamental  frequency  in  each  case  is  almost  com¬ 
pletely  unaffected.  The  higher  frequencies  are  affected  by  the  coupling  be¬ 
tween  flapwise  bending  and  torsion,  and  between  all  three  types  of  deformation 
when  twist  is  introduced. 

The  combined  effects  of  rotation  and  twist  on  blade  No.  1  are  illustrated 
in  Fig.  2.5o  In  the  untwisted,  nonrotating  case  the  fundamental  mode  is  iden¬ 
tified  as  predominantly  flapwise  bending;  the  second  mode  is  uncoupled  chord- 
wise  bending;  and  the  third  and  fourth  modes  are  coupled  flapwise  bending  and 
torsion.  The  effect  of  rotation  and  twist  is  to  couple  the  first  two  modes. 

A  comparison  of  the  results  in  Fig.  2. 5  ^ith  results  in  Ref.  1  for  a  beam  with 
the  same  bending  properties  but  with  torsion  neglected  shows  that  essentially 
no  change  has  been  introduced  by  the  presence  of  torsion.  The  fourth  frequency 
in  Fig.  2.5  differs  slightly  from  the  third  frequency  for  the  beam  in  Ref.  1, 
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TABLE  I 


BEAM  PROPERTIES 


Beam  No. _ 1 _ 2 


i — 

h 

— ► 

T 

— -  L  - 

h/L 

2.285-10-1 

6.225-10-2 

h/R 

l.OOO-lO-2 

l.OOO-lO'2 

e 

0.25  L 

0.15  L 

eo 

0 

0 

L 

0.015522 

0.04926 

pip 

0 

0 

72 

1.000 -10-1 

1.000 -10-2 

Ii[in.4] 

2.355  h3t 

8.1987  h3t 

I2[in.4] 

(2.355-10)h3t 

(8.1987-102)h3t 

Bx[in.6] 

32.692  h5t 

15.345  h5t 

B2[in.5] 

0 

0 

GJ| 

1.153 

1.405 

h 

2.190-10-5 

2.402-10'5 

h 

5-832-10-3 

7.114-10"3 

*Aluminum  is 

assumed. 
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while  the  third  frequency  in  Fig.  2.5  is  a  new  one  introduced  by  the  presence 
of  torsion.  It  is  seen  that  twist  has  little  effect  on  the  third  and  fourth 
frequencies  shown  in  Fig*  2 ,5 .  ,A  small  difference  between  the  results  in  Kef. 

1  and  the  present  results  is  introduced  by  the  inclusion  of  rotary  inertia 
in  the  present  analysis . 

The  results  for  blade  No.  2  presented  in  Fig  2,6  show  that  the  funda¬ 
mental  frequency  is  essentially  unchanged  by  the  presence  of  torsion  when  com¬ 
pared  with  results  for  a  similar  blade  reported  in  Ref.  h  For  this  blade  the 
fundamental  mode  is  predominantly  flapwise  bending..  The  second  and  third  modes 
for  the  untwisted,  nonrotating  blade  are  coupled  flapwise  bending  and  torsion, 
and  the  fourth  mode  is  uncoupled  chordwi.se  bending.  When  rotation  and  twist 
are  added,  the  three  higher  modes  exhibit  considerable  coupling,  and  it  be¬ 
comes  difficult  to  reach  any  general  conclusions.  When  compared  with  the  re¬ 
sults  in  Ref,  1  for  a  similar  blade  without  torsion,  it  is  seen  that  the'  ef¬ 
fect  of  the  presence  of  torsion  is  to  introduce  a  new  frequency  and  to  modify 
the  other  two  frequencies  a  moderate  amount.  That  these  two  frequencies  are 
not  modified  more  by  the  presence  of  torsion  is  to  be  expected  since  the  coupl¬ 
ing  for  this  blade  (and  also  for  blade  No,  l),  as  represented  by  the  amount 
of  offset  between  the  mass  and  elastic  axes,  is  relatively  small. 

The  neglect  of  rotary  inertia  has  a  negligible  effect  on  the  natural  fre¬ 
quencies  except  for  those  cases  in  which  there  is  chordwise  bending.  For  exam¬ 
ple,  in  the  case  of  the  untwisted  beam  No,  2  the  only  frequency  which  is  ap¬ 
preciably  affected  is  the  uncoupled  chordwise  bending  frequency,.  The .  magnitude 
of  this  effect  is  shown  in  Fig,  2,6, 

The  effects  of  centrifugal  force  coupling  on  beam  No.  2  are  shown  in 
Fig,  2.7«  Curves  with  and  without  centrifugal  force  coupling  are  shown  for 
©  =  0°  and  50° •  Only  the  9=0°  case  is  shown  for  the  second  frequency  to 
avoid  confusion  in  plotting.  The  ©  =  30°  case  for  the  second  frequency  is 
modified  by  a  slightly  smaller  amount.  The  9=0°  case  which  represents 
uncoupled  chordwise  bending  (the  fourth  frequency  for  the  nonrotating  beam) 
is  essentially  unaffected,  as  is  the  ©  =  30°  case  for  the  third  coupled  fre¬ 
quency,  It  can  be  seen  from  these  results  that  centrifugal  force  coupling 
can  have  an  appreciable  effect  on  some  of  the  vibration  characteristics. 
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3.  SIMPLE  MODEL  ANALYSIS 


SOBOLS 
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i0 

If 
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Mp 

Mp 
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r 

-> 

R 

t 

U 

-> 

v 

v 


functions  defined  immediately  following  Eq.  (3.13) 

offset  of  mass  c.g.  from  supporting  rod,  positive  forward 

nondimens ional  form  of  e,  e/r 

unit  vectors  along  the  x,y,z  axes  respectively 

moment  of  inertia  of  mass  m  about  supporting  rod 

moment  of  inertia  of  mass  m  about  its  own  c.g. 

moment  of  inertia  of  flywheel 

kinetic  energy  of  system 

stiffness  of  bending  spring 

stiffness  of  torsion  spring 

mass 

shaft  torque 

nondimens  ional  form  of  Mp,  Mp/mr2^2 
differential  operator 

length  of  supporting  rod  from  shaft  to  mass  m 
radius  vector  from  origin  to  element  of  mass  dm 
time 

potential  energy  of  system 
velocity  vector  of  mass  element  dm 
magnitude  of  v 
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axes  fixed  to  supporting  rod  and  mass  assembly 
y-coordinate  of  mass  element  dm 
stationary  axes 

orientation  angle  of  bending  hinge  axis 
angle  simulating  built-in  twist 
built-in  coning  angle 

phase  lag  of  motion  in  ^-coordinate  relative  to  motion  in  0- 
coordinate 

elastic  displacement  about  bending  hinge 

pseudo -static  value  of  0 

departure  of  0  from  0S 

amplitude  of  0,  also  initial  value  of  © 

nondime n s i onal  radius  of  gyration  of  mass  m  about  supporting  rod, 
\T E/mr2 

nondimens i onal  radius  of  gyration  of  mass  m  about  its  own  c.g., 
nT L0/mr2 

nondimensional  form  of  t,  Ht 

elastic  displacement  about  torsion  hinge 

pseudo -static  value  of  l 

departure  of  f  from 

initial  value  of  <( 

initial  value  of  <j' 

amplitude  of  j 

angular  displacement  of  shaft 
angular  velocity  vector  of  x,y, z  frame 
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components  of  oi  along  the  x,y,  z  axes,  respectively 
go  natural  frequency  of  characteristic  oscillation 

a>i,ac£  first  and  second  natural  frequencies  of  characteristic  oscillation 

Cl  rotational  velocity  of  shaft 

DESCRIPTION  OF  THE  MODEL 

In  order  to  examine  some  effects  of  nonlinearity  and  Coriolis  forces  in 
the  free  vibrations  of  a  rotating  elastic  blade  in  coupled  bending  and  torsion 
and  to  consider  the  effects  of  certain  parameters  on  the  static  deformation  of 
the  rotating  blade,  a  simple  model  with  a  small  number  of  degrees  of  freedom 
is  set  up  and  analyzed. 

The  model  consists  of  a  rigid  weightless  rod  on  one  end  of  which  is  mounted 
a  mass  and  the  other  end  of  which  is  connected  to  a  rotating  shaft.  The  connec¬ 
tion  to  the  shaft  is  through  a  hinge  with  axis  normal  to  the  rod  and  set  at  an 
angle  to  the  shaft.  A  spring,  restraining  motion  about  this  hinge,  simulates 
bending  stiffness.  In  addition,  the  rod  is  free  to  rotate  about  its  own  axis 
against  the  action  of  a  spring,  which  simulates  torsional  stiffness.  The  mass 
is  assumed  to  be  distributed  along  a  line  normal  to  the  rod,  simulating  the 
major  principal  axis  of  a  blade  cross  section,  with  its  center  of  gravity  dis¬ 
placed  from  the  rod,  simulating  an  offset  of  the  mass  axis  of  the  blade  from 
the  elastic  axis . 

The  orientation  of  the  model  relative  to  a  set  of  fixed  axes  and  the  gen¬ 
eralized  coordinates  defining  its  configuration  are  shown  in  Fig.  3«1*  The 
final  orientation  is  reached  by  aligning  the  model  initially  with  the  fixed 
axes  and  then  executing  a  sequence  of  rotations.  The  fixed  axes  x^,,  y^,  z^ 
form  an  orthogonal  set  oriented  so  that  the  xp-axis  is  coincident  with  the 
shaft  centerline .  Their  origin  is  at  the  intersection  of  the  rod  and  the 
shaft  centerline  and  is  coincident  with  the  origin  of  the  model  axes  x,  y,  z. 

The  x-axis  lies  along  the  rod,  the  y-axis  is  parallel  to  the  line  along  which 
the  mass  lies,  and  the  z-axis  completes  the  orthogonal  set. 

The  model  is  initially  aligned  so  that  the  x,  y,  z  axes  are  coincident 
with  the  xp,  yp,  zp  axes,  respectively.  The  following  rotations,  positive 
in  the  right-handed  sense,  of  the  x,  y,  z  frame  are  then  executed  in  sequence: 

1.  A  rotation  about  the  Zp-axis  through  the  angle  \|r  to  the  position  xi, 
yi,  zi.  \|r  then  defines  the  shaft  rotation. 

2.  A  rotation  about  the  yi-axis  through  the  angle  -7  to  the  position  X2, 

y2,  z2 •  -7  then  defines  a  built-in  coning  angle. 
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Bending  Hinge  Axis 


3.  A  rotation  about  the  x2-axis  through  the  angle  a  to  the  position  x3, 
y3,  z3 •  The  y3-axis  then  defines  the  position  of  the  hinge  axis. 

4 .  A  rotation  about  the  y3-axis  through  the  angle  -0  to  the  position  X4 
y4,  Z4 .  This  represents  a  rotation  about  the  hinge  axis  simulating 
bending  displacement.- 

5.  A  rotation  about  the  X4~axis  through  the  angles  p  and  <j>  in  sequence 
to  the  final  position  x,  y,  z.  The  angle  f3  simulates  built-in  twist 
and  the  angle  <f>  elastic  twist. 

The  angles  a,  P,  and  7  are  constants  and  constitute  parameters  in  the 
problem.  The  angles  \|r ,  ©,  and  are  generalized  coordinates  representing 
shaft  rotation,  bending,  and  torsional  displacement,  respectively. 

DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

The  equations  governing  the  motion  of  the  model  are  now  derived  using 


Lagrange's  equation.  Toward  this  end  it  is  necessary  to  obtain  an  expression 
for  the  kinetic  energy  of  the  system  in  terms  of  the  generalized  coordinates. 

Assuming  a  flywheel  of  moment  of  inertia  Ip  to  be  mounted  on  the  shaft; 
and  defining  m  as  the  magnitude  of  the  mass  mounted  on  the  rod,  the  kinetic 
energy  of  the  system  may  be  written, 


K  = 


(3-1) 


where  v  is  the  magnitude  of  the  velocity  vector  v  of  an  element  of  the  mass  m. 
v  may  he  developed  from  the  relation, 


■+  -4  -4 

V  =  03  X  R  (3.2) 


where  cu  is  the  angular  velocity  vector  of  the  x,  y,  z  frame  and  $  is  the  radius 
vector  of  dm.  Substituting 


®  =  “X1  +  “yJ  +  (3-3) 

R  =  ri  +  yj  (3.4) 

where  i,  j,  k  are  unit  vectors  along  the  x,  y,  z  axes,  respectively,  into  Eq. 
(3.2)  the  following  is  obtained, 

-4-4  -4 

v  =  -  y^z1  +  rU3zJ  +  (yubc-r<Uy)k  •  (3.5) 

Thus, 


v‘  =  +  (yu^x-rcDy) 


2.2 


(3-6) 


and  Eq.  (3*1)  may  now  be  written, 


K  ~  2  +  I  +  2  I(ubc+<«z)  -  mreo)x  u>y  (3-7) 
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Neglecting  gravity  forces,  the  potential  energy  of  the  system  may  be  written 
as  follows, 


=  \  k0®2  +  |  ks/ 


where  kQ  and  are  the  spring  constants  of  the  springs  restraining  motion  in 
0  and  </>  coordinates . 

Substitution  of  Eqs .  (3*7);  (3-8)  and  (3-9)  into  Lagrange's  equation, 


i.  (*)  .  j*  .  .  o  (1=1,2, 5) 

dt  \dqj_  /  dq*  dq* 


(3.10) 


where  qi  =  0,  q2  =  ‘t>  03  =  yields  the  following  differential  equations. 


[mr2+I  sin2 


(p+^)]0  +  [rare  cos(p+^)]^  +  -mr2  sin  a  cos  y+mre  a2  cos(p+^) 


+  I  sin(p+jz(){  -sin(|3+^)sin  OL  cos  7+ ai  cos(p+^))J  ^ 

+  |~2 1  sin(p+jzO  {cos(p+5z( )sin  a  cos  7+  ai  sin(p+jzO  )-2mrea2  sin(f3+^)J  ^ 
+  [21  sin(p+^)cos(g+^)  ]Q</>  +  [-rme  sin((3+sz0  ]p 
+  jmr2aia2+mre{ai  cos(p+^)sin  a  cos  7+(ai-af  )sin(  $+</> ) } 

-  lag  sin(p+^)  (cos(p+^)sin  a  cos  7+ai  sin(p+$zO}J 


+  kQ©  =  0 


(3.H) 


[mre  cos(p+^)]8  +  +  jla2-mre (cos(p+jOsin  a  cos  7+ ax  sin(p+^)}J  \jr 

+  |2mrea2  sin(p+^)+2I  sin(p+$^)  {cos(p+^)sin  (2  cos  7+ ai  sin(p+^))J  ©^ 

+  [-1  sin(p+^)cos(p+$z()  ]©2 
+  jmrea2 (-sin(p+^)sin  a  cos  7+ax  cos(pV)) 

+  I(-sin(p+$zf  )sin  Ot  cos  7+ai  cos  (p+^) }  (cos(p+0  )sin  a  cos  7+ai  sin(p+^)}J 


+  k  A  =  0 


(3.12) 


jj-mr2sin  a  cos  7+mrea2  cos(p+^)+I  sin( P+^ ) { -sin( fi+l )sin  Ot  cos  7+ai  cos  (p+^)}j© 
+  |la2  -mre[cos(p+^)sinacos  7+ai  sin(p+^))J  *</> 

+  jlF+mr2(sin2  ot  cos2  7+a2  )-2mrea2  [cos(p+^)sin  ot  cos  7+ai  sin(p+$zf)} 

+  Ia2+I{-sin(p+$z()sin  ot  cos  7+ai  cos(p+^)}2J  \jr 

+  |~2mr2aia2+2mre (a|  sin(p+^)-ai  cos(p+^)sin  a  cos  7-aia2  sin(p+^)} 

+  2Iaia2-2Ia2  cos(p+^) {-sin(p+^)sin  a  cos  7+ai  cos(p+  ^))J  ^ 

+  £-2rnrea2(-sin(p+^)sin  ot  cos  7+ai  cos(p+^)) 

-  21  (cos ( p+^ )sin  Ot  cos  7+ai  sin(p+^ ) ) ( -sin( p+  ^)sin  a  cos  7+ai  cos  ( p+^ ) )] 

+  [21  cos(p+szf)(  -sin(p+jz$)sin  a  cos  7+ax  cos(p+^)}J  ©$ 

+  [mre-ai  cos(p+^)-Ia2  sin(p+^)cos(p+^)  ]©2 


jmre {sin(p+^)sin  Ot  cos  7~ai  cos(p+^)}J  ft 


(3*13 


kl 


where 


ai  =  -sin  9  sin  7  +  cos  0  cos  Ot  cos  7 


a2  =  cos  9  sin  7  +  sin  9  cos  a  cos  7 


SPECIALIZATION  TO  THE  CASE  OF  CONSTANT  SHAFT  SPEED 

The  problem  is  now  specialized  to  the  case  of  constant  rotational  velocity 
of  the  shaft  by  setting  =  0,  if  =  ft,  and  the  equations  are  put  Into  a  nondimen- 
sional  form  by  defining  the  nondimens ional  parameters , 

P  = 


=b  = 

“t  - 

and  introducing  the  nondimensional  time  variable , 

t  =  fit 


It  is  seen  that  p  is  the  nondimensional  radius  of  gyration  of  the  mass  about 
the  rod  axis,  e*  is  the  nondimensional  offset  of  the  mass  center  of  gravity  from 
the  rod  axis  and  and  are  respectively  the  nondimensional  natural  fre¬ 
quencies  in  restrained  bending  and  restrained  torsion  when  the  shaft  is  not  ro¬ 
tating. 

Division  of  Eq.  (3.1l)  by  ^mr2  and  Eq.  (3. 12)  by  ft2I  now  yields, 
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[1+p  sin2(p+^)]@"  +  e  cos(p+^)^M 

+  j^p2  sin(pV){cos(pV)sin  a  cos  7+ax  sin(p+^) }-2ea2  sin(p+^) 
+  2p2  sin(P+jZ<)cos(p+^)9  V  -  e  sin(P+/)jzf'2 
+  laxas+efai  cos(p+j^)sin  a  cos  7+(a2-a|  )sin(p+^)  ] 

-  p2a2  sin(pV)  (cos(0+jOsin  (X  cos  7+ai  sin(0+^)}J 


+  o^Q  =  0 


(3-14) 


4  +  Zs  cos ( 0+^)0 "  +  [2  §5-  a2  sin(p+szO 

P  M 

+  2  sin(0+j/)  {cos(0+^)sin  a  cos  7+ai  sin(0+/))T  O' 

-  sin(p+/)cos(p+^)©,2+  [jp-  a2(-sin(0+jzOsin  a  cos  7+ai  cos(0+^)} 

+  (-sin(0+^)sin  a  cos  7+ai  cos(p+/) }  (cos(p+^)sin  a  cos  7+ai  sin(0+^))J 
+  aS.$  =  0 


(3-15) 


where  primes  denote  differentiation  with  respect  to  r,  and  a!  and  a2  are  as 
defined  in  the  preceding  section. 

Recognizing  that  constant  shaft  speed  represents  the  limiting  case  of  in¬ 
finite  flywheel  inertia,  the  term  Ipi|J  in  Eq.  (3.13)  can  he  seen  to  remain  finite 
and  equal  to  the  shaft  torque,  which  may  then,  from  Eq.  (3.I3),  he  written  in 
the  following  nondimens ional  form. 


%  =  £“Sin  OL  cos  7+ea2  cos(p+^) 

+  p2  sin(p+^)  (-sin(p+^)sin  &  cos  7+ai  cos(p+^)]J  ©" 

+  jp2a2-*e(cos(p+^)sin  OL  cos  7+ai  sin(p+^))J 

+  |“2aia2+2e  (a|  sin(p+^)-ai  cos(p+^)sin  a  cos  7-aia2  sin(p+{z()} 

+  2p2aia2*2p2a2  cos(p+^)  {-sin(p+^)sin  OL  cos  7+ai  cos  0+/))]  O' 

+  jj-2ea2(-sin(p+^)sin  Q!  cos  7+  ai  cos(p+?0} 

-  2p2{cos(P+^)sin  a  cos  7+ai  sin(p+pO)  (-sin(p+$0sin  a  cos  7+a!  cos(pV))]^ 
+  2 p2  eos(p+^)  [-sin(p+^)sin  a  cos  7+ai  cos(p+^)]©y' 

+  [eai  cos(3+jz()-p2a2  sin( p+^ )cos ( p+$zf )  ]© ,2 

+  e  [sin(p+^)sin  a  cos  7-ai  cos(p+^)]^’2  (3.16) 


where 


Mp 


Mrp 

mr2^2 


and  Mp  is  the  dimensional  torque. 

To  facilitate  solution,  it  is  desirable  to  rearrange  Eqs.  (3.14)  and  (3. 15) 
in  the  form, 


fi©"  +  fW"  =  -f3e-/'  -  f 4/'  -  f5iz(t2  -  5^0  -  fe  (5.17) 

f7©’’  +  ft'  =  -f80 '  -  f9©'2  -  33^  -  f1D  (5.18) 


where 


bb 


fi  =  1  +  p2  sin2(p+$zO 

f 2  =  .  e  cos(p+^) 

f3  =  2p  sin(p+^)cos(p+^) 

f4  =  -2p2  sin(0+fO {cos(0+^)sin  a  cos  7+ax  sin(p+^) )-2ea2  sin(0+fO 

f5  =  -e  sin(p+$zO 

f6  =  aia2  +  e[ai  cos (p+jzO sin  Ot  cos  y+( af-af  )sin(p+$zi) } 

-  p2a2  sin(p+^)  (cos(pV)sin  a  cos  7+ai  sin(0+^)) 

f7  =  %  cos(p+^) 

P 

f8  =  2  Jg-  a2  sin(P+^)+2  sin(p+^ )  {cos(  )sin  ot  cos  7+ai  sin(p+^)} 

P 

f g  =  -sin(p+^)cos(p+^) 

fi^  =  ro  a2 [-sin(0+^)sin  a  cos  7+ax  cos(f3+jz0) 

P2 

+  (-sin(p+^)sin  a  cos  7+ax  cos(P+jO  }(cos(0+^)si'n  Ot  cos  7+ax  sin(0+$O) 

Solving  Eqs<>  (3-17)  and  (3-l8)  for  0"  and  <$"  in  terms  of  0  and  f  and  their 
first  derivatives  yields  the  differential  equations  in  the  following  form, 

0"  =  - - -  [-faO'^'-f^'-M^-^e-fR+fafse’+fafse^+fa^+fafxo]  (3-19) 

fi-f2f7  b 

</>"  =  [-f  if 8e  ’  -f if©6 ,2-f  icE|/-f  xf  10+f Tfis©  ’  / ' 

+  f7f4^,+f7fs^’2+f7^©+f7f6J  •  (3-20) 

Equations  (3°19)  and  (3*20)  are  now  in  suitable  form  for  solution  on  a 
digital  or  analog  computer. 
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SOLUTION  OF  THE  PSEUDO-STATIC  PROBLEM 


It  is  of  interest  to  determine  the  static  configuration  of  the  rotating 
model,  that  is,  the  static  displacements  under  the  action  of  centrifugal 
forces.  This  problem  may  be  termed  the  pseudo-static  problem.  Its  solution 
permits  the  setting  up  and  solution  of  linearized  differential  equations  for 
small  motions  about  the  pseudo-static  configuration. 

The  appropriate  equations  are  obtained  by  eliminating  all  terms  contain¬ 
ing  derivatives  of  0  and  jzi  from  Eqs.  (3-17)  and  (3.18),  yielding, 


<u§es  +  f6  =  0 

(5-21) 

+  ^10  =  0  . 

(3.22) 

These  equations  are  nonlinear,  with  f6  and  fiG  being  transcendental  func¬ 
tions  of  the  dependent  variables.  Since  it  is  not  feasible  to  obtain  an  ana¬ 
lytical  solution  in  closed  form,  the  following  iterative  procedure  was  applied 

Equations  (3-2l)  and  (3-22)  are  linearized  with  respect  to  departures  A© 
and  A^  from  trial  values  9n  and  ^n,  respectively,  of  the  variables,  yielding, 

^0n  +  ^AO  +  f6n  +  (|^)  ^  AO  +  (|^)n  a/  =  0  (3.23) 

KK  +  ^  *  (^)  n  *  '  0  (5.21,) 

where  subscript  n  denotes  values  at  9  =  ©n,  </>  =  ^n. 

Equations  (3*23)  and  (3-24)  are  now  rearranged  in  the  form, 


+  fl2n^  = 

-f6n  -  ^©n 

(5-25) 

fl3n^®  +  = 

-fion  - 

(3-26) 

where 
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fi: 


_  ^f  6 

S9 

=  af  -  a§  -  e(a2  cos (p+^) sin  a  cos  7+4aia2  sin  (p+jzO} 
-  p2  sin(p+^)(a*i  cos(P+^)sin  a  cos  7+(af-a|  )sin(p+jzO  } 


fi*  = 


2  2 


ef-ax  sin(pV)sin  Of  cos  7+(ax-a2  )cos(p+$zO } 

p2a2  £[cos2(p+^)-sin2(P+^)  }sin  a  cos  7+2ai  sin(  p+^)cos  ( p+^  )J 


df 


13 


IQ. 


a© 


f*14 


=  —  (-ai  sin(p+^)sin  a  cos  7+(ai-a2 )cos(p+^ ) } 

P 

-  a2  J^{cos2(p+^)-sin2(p+^)  ]sin  a  cos  7+2ai  sin( pV)cos  ( P+^)J 
_  10 

-  -  a2 (cos (p+^)sin  a  cos  7+ai  sin  (p+^)} 

(cos(p+^)sin  a  cos  7+ai  sin(p+^))2 

+  [-sin(p+^)sin  a  cos  7+  ax  cos(p+^)}2 

Solution  of  Eqs *  (5-25)  and  (5-26)  for  A©  and  yields. 


A© 


-e4ne2n+e5nfl2n 


e3. 


A^  =  "esnein+e4iifi3n 


e3 


n 


where 


(3-27) 

(3-28) 
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el 

=  K  + 

fll 

e2 

+ 

ii 

fl4 

e3 

=  ele2 

-  fl2fl3 

e4 

■ 

+  f  6 

e5 

II 

+  fio 

Equations  (3*27)  and  (3*28)  may  "be  applied  in  conjunction  with  the  iter¬ 
ation  formulae , 


en+l 

=  ©n  +  A© 

(3-29) 

^n+1 

=  +  a/ 

(3-30) 

using  as  initial  values,  ©i  =  0,  =  0.  The  process  has  been  found  to  con¬ 

verge  rapidly  in  the  cases  that  have  been  considered  in  the  present  work. 


FORMULATION  AND  SOLUTION  OF  THE  LINEARIZED  EQUATIONS 

In  order  to  assess  the  significance  of  nonlinear  effects  in  the  problem 
under  consideration,  it  is  desirable  to  obtain  also  solutions  to  linearized 
equations  for  small  perturbations  ©  and  jzf  from  the  pseudo -static  configuration 

Application  of  small  perturbation  theory  to  Eqs.  (3-17)  and  (3.18)  yields 


fise"  +  (5|+  fiis)e  +  f ZJ"  +  Uj'  +  fizj  =  o  (3.31) 


f7sQ"  +  f8s©’  +  fi3Se  +  i"  +  (a^+fi4s)^  =  0 


(3-52) 


in  which  subscript  s  denotes  values  corresponding  to  the  pseudo-static  config¬ 
uration  ©  =  ©g,  </>  =  jQ. 

Putting  Eqs.  (3»3l)  and  (3-32)  into  operator  form,  using  symbol  p  to  de¬ 
note  the  differential  operator,  and  expanding  the  determinant  of  coefficients, 


48 


the  following  characteristic  equation  is  obtained, 


Cip4  +  c2p2  +  c3  =  o  (3-33) 

where 

Cl  =  fls  -  f2sf7s 

C2  =  fls(u^+fl4s)  +  0%  +  flls  -  f2sfl3s  -  f7sfl2s  "  f4sfss 
C3  =  (c^+flls)(<*^+fl4s)  "  fl2gfl3s 

The  terms  in  p  and  p3  are  seen  to  vanish. 

The  roots  of  this  equation  are 

pf  =  gi-  (-c2  +  >/ci-4cic3) 

pi  =  2^-  (-c2  -  '/c|-4c1c3) 

and  the  characteristic  frequencies  are  given  by 


0)i 


(3.34) 


(3.35) 


The  characteristic  mode  shapes  may  be  determined  by  assuming  a  solution 
of  the  form, 


©  = 


©0  COS  0)T 


(3-36) 


<f>  =  <j)Q  COS 


CDT 


J  t 

+  gp  sin  cdt 


=  i 


1  cos 


(CDT-S  ) 
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(3-37) 


introducing  Eqs.  (306)  and  (3*57)  into  Eq.  (3*31 )>  and  equating  the  sum  of 
the  coefficients  of  the  cos  cot  and  the  sin  cot  terms  respectively  to  zero,  the 
following  result  is  obtained, 


4o/®0  = 

(  f2s032 -f  12s  )  ( -f  Igof-td^+f  ns  ) 

(3-38) 

(f2saf-f12s)2  +  f|s0? 

4'of^o 

f4saj2  ( -f  istiP-Kl^+f  ns ) 

(f2suf  -fi2s)2  +  f4s® 

(3-39) 

where  co  =  cox,  o^. 

The  mode  shapes  may  he  expressed  alternatively  in  the  form, 


(3.40) 

(3-41) 


where,  from  Eq..  (3*37)  it  is  seen  that  4x /©0  is  the  relative  amplitude  of  dis- 
placements  in  the  two  coordinates  and  e  is  the  phase  lag  of  the  oscillation  in 
the  ^-coordinate  relative  to  that  in  the  ©-coordinate. 

A  solution  involving  only  one  characteristic  mode  of  oscillation  may  he 
obtained  by  selecting  as  initial  conditions, 


©  =  es  +  oG  ,  4  =  4S  +  40  ,  </>'  =  4'0  (3-42) 


where  ©0  may  be  selected  arbitrarily  within  the  limitations  imposed  by  the  as¬ 
sumption  of  small  perturbations,  and  $z40  and  fa  are  then  determined  from  Eqs. 
(3.38)  and  (3-39)- 

It  should  be  noted  that  the  existence  of  a  phase  difference  between  oscil 
lation  in  the  two  coordinates  is  associated  with  the  presence  of  the  terms 
f 4^*  and  f qs©^  in  Eqs.  (3-31)  and  (3-32)  respectively.  These  terms  originate 
in  the  terms  in  6^  and  in  Eqs.  (3*11)  and  (3-12),  which  are  due  to  the  pres 
ence  of  Coriolis  forces. 


DISCUSSION  OF  RESULTS 

A  series  of  computations  on  the  simple  model  were  performed  using  an  auto 
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matic  digital  computer.  There  computations  were  limited  to  the  case  of  con¬ 
stant  rotational  velocity  of  the  shaft  and  zero  "built-in  coning  angle  (7  =  0). 

The  pseudo -static  configuration  was  determined  hy  means  of  the  iterative 
procedure  developed  earlier ,  and  corresponding  characteristics  of  the  linear¬ 
ized  system  for  small  perturbations  from  this  configuration  were  computed.  In 
each  case  additional  computations  were  performed  in  which  the  terms  f4S^'  and 
f8s0’  in  Eqs.  (3*31)  and  (3*32)  were  omitted.  As  discussed  previously,  these 
te'rms  represent  the  influence  of  Coriolis  forces,  so  that  a  comparison  of  re¬ 
sults  obtained  with  and  without  their  inclusion  provides  a  means  of  assessing 
the  importance  of  the  Coriolis  forces. 

These  results  are  presented  in  Figs.  3*2  to  3*8  inclusive.  Figures  3.2 
and  3.3  show  the  effect  of  varying  the  mass  offset  parameter  e,  with  the  param¬ 
eter  p0,  representing  the  nondimens ional  radius  of  gyration  of  the  mass  about 
its  center  of  gravity,  and  the  parameters,  a,  j3,  and  being  maintained 
constant.  Since  the  parameter  p  must  be  varied  accordingly,  the  maintenance 
of  a  constant  value  for  implies  that  the  variation  of  e  does  not  involve 
merely  a  shifting  of  the  mass  relative  to  the  supporting  arm  but  involves  also 
changes  in  m  or  or  both.  The  value  of  p0  selected  for  this  case  represents 
a  rather  extreme  value,  applicable  to  a  short,  wide  blade. 

Figure  3.2  shows  the  substantial  pseudo-static  deformation  occurring  in 


Fig.  3.2.  Effect  of  mass  offset  on  pseudo-static 
displacements  of  model.  a  =  30°,  p  =  15 ° ,  u>b  = 

1/3,  =  1,  Po  =  0.1732. 
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this  case.  It  should  he  noted  that  the  static  twist  decreases  with  increase 
in  offset^  when  the  center  of  gravity  of  the  mass  is  behind  the  elastic  axis. 
This  occurs  despite  the  fact  that  the  relative  values  of  the  moment  of  inertia 
and  torsional  stiffness  about  the  supporting  arm  remain  the  same  because  of 
the  constancy  of  a>t,  which  fact  implies  that  centrifugal  twisting  moment,  be¬ 
fore  deformation,  remains  the  same.  It  must  be  concluded  that  the  variation 
in  twist  is  associated  with  a  component  of  centrifugal  force  normal  to  the 
coning  surface  on  which  the  supporting  arm  revolves.  This  effect  is  intro¬ 
duced  through  the  term  in  mre^2  in  Eq.  (3*12)  and  terms  deriving  from  it  in 
later  forms.  It  has  been  called  "centrifugal  force  coupling"  in  Ref.  2,  and 
shown  there  to  have  a  substantial  effect  on  natural  coupled  frequencies  of  vi¬ 
bration.  In  the  present  case,  since  7  =  0,  the  coning  of  the  supporting  arm 
is  associated  solely  with  the  displacement  ©s.  With  positive  ©s  and  a  positive 
value  for  (a+| 3+$zfs),  this  effect  opposes  that  of  centrifugal  twisting  moment. 

It  can  be  expected  to  be  more  pronounced  in  the  case  of  blades  with  built-in 
coning  angle. 

Figure  3*3  shows  the  effect  of  mass  offset  on  the  natural  vibration  char¬ 
acteristics  of  the  system  linearized  with  respect  to  the  pseudo-static  config¬ 
uration.  As  can  be  expected,  it  is  seen  that  the  increased  coupling  between 
bending  and  torsion  associated  with  increasing  mass  offset  separates  the  natural 
frequencies  and  alters  the  natural  mode  shapes. 

It  is  seen  also  that  the  Coriolis  forces  introduce  substantial  phase  dif¬ 
ferences  between  motion  in  the  two  coordinates,  particularly  in  the  case  of 
the  first  or  predominantly  bending  mode,  where  the  phase  angle  is  large  through¬ 
out  the  range  of  e  considered.  In  the  case  of  the  second  mode,  where  torsional 
motion  predominates,  the  phase  angle  is  substantial  only  at  small  values  of  e. 
When  e  is  zero  the  only  coupling  between  bending  and  torsion  is  through  the 
Coriolis  forces,  and  the  phase  difference  is  then  90%  /  leading  ©  by  this 
amount  in  the  case  of  the  first  mode  and  lagging  by  this  amount  in  the  case  of 
the  second  mode.  Furthermore,  the  Coriolis  forces  are  seen  to  have  a  substan¬ 
tial  effect  on  the  mode  shape  of  the  first  mode  and  a  somewhat  modest  effect 
on  the  corresponding  frequency.  The  corresponding  effects  on  the  second  mode 
and  frequency  are  seen  to  be  small  or  negligible.  It  should  be  noted  here  that 
the  apparent  absence  in  some  cases  of  curves  associated  with  neglect  of  Coriolis 
forces  is  explained  by  the  fact  that  such  curves  are  indistinguishable  from 
the  corresponding  solid-line  curves,  and  the  effect  of  these  forces  is  thus  very 
small . 

Figures  3.4  and  3.5  provide  results  corresponding  to  those  of  Figs.  3.2 
and  3*3  f*or  a  different  case,  namely  one  involving  a  much  smaller  value  of 
p0  and  consequently  more  realistic  in  relation  to  propeller  or  helicoptor 
rotor  blades.  Similar  trends  are  observed,  except  that  Coriolis  force  effects 
are  considerably  reduced,  but  still  substantial  with  respect  to  phase  differ¬ 
ences  in  the  first  mode. 
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Fig.  3.3.  Effect  of  mass  offset  on  natural  vibration  characteristics 
of  model,  a  =  3O0,  p  =  -15°,  0%  =  1/3,  *t  -  1,  p0  «  0.1732. 


Fig.  3-4.  Effect  of  mass  offset  on  pseudo-static 
displacements  of  model,  a  =  3O0,  p  =  -15%  = 

1/3;  a>t  =  1,  p0  =  0.05. 

Figure  3*6  shows  the  effect  of  varying  the  tending  hinge  orientation 
angle  while  maintaining  the  orientation  of  the  principal  axis  of  the  mass  fixed. 
This  involves  varying  a  and  p  so  that  a+p  remains  constant ,  and  simulates  a  sit¬ 
uation  in  which  mean  "blade  angle  is  kept  constant  while  built-in  twist  is  varied. 
All  other  parameters  were  maintained  constant.  Curves  of  ©s  and  are  not 
shown,  as  variations  in  those  parameters  were  small.  For  a  variation  of  a  from 
15 °  to  45°,  ©s  varied  from  O.780  to  0-97°  and  varied  from  -6.5O0  to  -6.6O0. 

It  is  seen  from  Fig.  3*6  that  first  mode  characteristics  are  affected  very  sub¬ 
stantially  by  changes  in  a,  the  phase  difference  between  the  ©  and  {  motions 
especially  varying  over  a  very  wide  range.  The  effect  on  second  mode  charac¬ 
teristics  is  much  smaller,  although  still  considerable. 

Figures  3*7  and  3*8  show  the  effect  of  varying  the  rotational  velocity  of 
the  shaft  while  other  parameters  remain  constant.  The  information  in  Fig.  3.7 
is  principally  of  value  in  estimating  the  pseudo -static  torsional  deformation 
corresponding  to  a  given  rotational  velocity.  This  deformation  can  be  expected 
to  depend  primarily  on  the  parameter  5^  in  the  case  of  a  blade  without  built- 
in  coning  angle,  although  from  results  discussed  earlier  it  can  be  seen  also  to 
depend  somewhat  on  the  parameters  e  and  25^.  From  Fig.  3.7  it  can  be  seen  that 
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PHASE  ANGLE,  €  DEGREES  PHASE  ANGLE,  -<=  DEGREES 


"b 

Fig-  3*7-  Effect  of  rotational  velocity  on  pseudo -static 
displacements  of  model.  o\/(X^  =  3,  a  =  45%  (3  =  “15°;  P  = 

0.2,  e  =  -0.1. 

the  torsional  displacement  will  exceed  20$  of  the  initial  blade  angle  (a+0)  if 
is  less  than  about  2,  that  is,  if  the  rotational  velocity  is  greater  than 
about  one-half  the  value  of  restrained  torsional  frequency  corresponding  to  zero 
rotational  velocity. 

Figure  3*8  indicates  an  increasing  prominence  of  torsion  relative  to  bend¬ 
ing  in  both  modes  as  rotational  velocity  becomes  large.  It  indicates  further 
a  marked  sensitivity  of  the  phase  difference  between  coordinates  in  the  first 
mode  to  variation  in  rotational  velocity,  at  least  in  a  limited  range  of  rota¬ 
tional  velocity.  The  phase  angle  is  seen  to  approach  zero  at  large  values  of 
rotational  velocity.  A  somewhat  different  situation  is  seen  to  exist  in  the 
case  of  the  second  mode,  where  the  phase  angle  increases  with  increase  in  ro¬ 
tational  velocity. 

Digital  computer  solutions  to  the  nonlinear  differential  equations  were 
also  obtained,  using  a  Runge-Kutta  procedure.  Initial  conditions  were  estab¬ 
lished  on  the  basis  of  natural  vibration  characteristics  determined  from  the 
linearized  equations,  that  is,  by  applying  Eqs.  (3-42),  using  values  from  Eqs. 
(3*38)  and  (3-39)  and  the  pseudo-static  displacements.  With  such  initial  con- 
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ditions,  the  linearized  system  responds  in  only  one  of  the  natural  modes ,  and 
comparison  with  the  corresponding  response  of  the  nonlinear  system  provides  a 
means  of  assessing  the  extent,  to  which  nonlinear  effects  distort  the  motion. 

Results  were  obtained  for  only  one  case  and  are  shown  in  Fig.  3*  9-  Figure 
3.9(a)  illustrates  the  response  of  the  first  mode  when  the  initial  bending  dis¬ 
placement  from  the  pseudo-static  configuration  is  10° .  Displacement  in  the 
bending  coordinate  predominates  in  that  mode  and  its  time  history  is  seen  to  be 
distorted  only  slightly  by  nonlinear  effects.  There  is  a  slight  increase  in 
period  and  a  very  slight  but  irregular  variation  in  amplitude.  The  torsional 
response  is  seen  to  be  strongly  influenced  by  nonlinear  effects  in  a  manner 
which  suggests  that  there  is  substantial  coupling  with  the  second  mode.  The 
slight  variation  in  amplitude  of  the  bending  motion  is  likely  associated  with 
this  coupling. 

The  response  in  the  second  mode  for  an  initial  bending  displacement  of  2° 
from  the  pseudo -static  configuration  is  shown  in  Fig.  3 • 9(t> ) -  In  this  case, 
displacement  in  the  torsional  coordinate  predominates  and  has  an  amplitude  of 
about  22.5°.  It  is  seen  that  there  is  jan  appreciable  increase  in  period  caused 
by  nonlinear  effects,  but  otherwise  only  a  slight  distortion  of  the  motion  in 
both  coordinates.  Again,  it  is  likely  that  this  distortion  is  due  to  coupling 
with  the  first  mode. 

The  solutions  were  not  carried  far  enough  to  ascertain  whether  there  is 
a  decay  or  divergence  of  the  oscillations.  The  fact  that  such  may  exist  is  not 
inconceivable,  in  view  of  the  fact  that  the  system  is  not  necessarily  conserva¬ 
tive.  It  has  been  seen  to  be  conservative  when  linearized  with  respect  to 
small  perturbations  from  the  pseudo-static  configuration.  However,  with  impo¬ 
sition  of  the  condition  of  constant  shaft  rotational  velocity  it  is  a  driven 
system,  and  it  is  possible  that  nonlinear  effects  may  result  in  a  transfer  of 
energy  to  or  from  it  through  the  shaft. 

The  results  obtained  indicate  that,  at  least  for  the  case  considered,  any 
such  divergence  or  decay  will  be  small  and  probably  represent  a  negligible  ef¬ 
fect  in  comparison  with  aerodynamic  effects  in  the  case  of  an  actual  blade. 

It  is  possible  that  a  different  choice  of  parameters  or  the  introduction  of 
built-in  coning  may  produce  a  different  result.  This  requires  further  inves¬ 
tigation. 

On  the  basis  of  the  present  results  it  appears  that  the  effect  of  Coriolis 
forces  is  likely  to  have  a  greater  practical  significance  than  the  effect  of 
nonlinearity,  particularly  since  it  does  not  depend  upon  the  existance  of  large 
motions.  This  relates  mainly  to  the  problem  of  blade  flutter,  since  the  flutter 
phenomenon  is  highly  sensitive  to  phase  differences  between  motion  in  bending 
and  torsion.  The  phase  differences  associated  with  the  presence  of  Coriolis 
forces  may  conceivably  alter  the  balance  in  the  flutter  problem  sufficiently  to 
change  the  conditions  for  flutter  significantly. 
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Solution  to  Nonlinear  Equations 
Solution  to  Linearized  Equations 


6o 


Solution  to  Nonlinear  Equations 
Solution  to  Linearized  Equations 


Other  effects  which  have  not  been  considered  in  the  present  study,  but 
which  can  be  expected  to  be  of  considerable  importance  in  some  cases,  are 
those  of  nonlinearity  of  the  torsional  spring  and  of  centrifugal  tension  on 
torsional  stiffness  and  on  pseudo-static  deformation,  as  discussed  in  the 
Introduction.  Their  introduction  into  the  present  analysis  should  not  re¬ 
sult  in  undue  complication  and  would  represent  an  appropriate  and  desirable 
extension  of  the  present  work. 
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k .  CONCLUDING  REMARKS 


A  practical  numerical  method ,  suitable  for  implementation  on  an  automatic 
digital  computer,  has  been  developed  for  determining  the  natural  vibration 
characteristics  of  twisted  rotating  and  nonrotating  blades  in  coupled  bend¬ 
ing  and  torsion.  A  limited  numerical  study  indicates  that  the  method  is  an 
efficient  one  for  including  the  effects  of  bending-torsion  coupling  and  pre- 
twist.  The  nature  of  the  coupling  is  complicated  and  a  much  more  extensive 
parametric  study  would  be  needed  in  order  to  draw  general  conclusions.  It 
can  be  said,  however,  that  centrifugal  force  coupling  can  have  an  appreciable 
effect  when  there  is  a  substantial  offset  of  the  mass  axis  from  the  elastic 
axis. 


In  order  to  investigate  some  effects  of  nonlinearity  and  Coriolis  forces 
in  the  rotating  blade  vibration  problem,  a  study  has  been  made  of  a  simple 
model  with  a  small  number  of  degrees  of  freedom.  Computations  performed  on 
this  model  indicate  the  following: 

(1)  Th6re  is  an  effect  of  centrifugal  force,  apart  from  the  familiar 
centrifugal  twisting  moment,  on  the  torsional  deformation  when  the  mass  axis 
of  the  blade  is  offset  from  the  elastic  axis.  It  may,  in  some  cases,  modify 
the  static  deformation  of  the  rotating  blade  substantially,  and  tends  to  in¬ 
troduce  additional  coupling  between  bending  and  torsion  when  the  blade  is 
vibrating,  as  discussed  also  in  the  case  of  the  continuous  blade. 

(2)  The  presence  of  Coriolis  forces  causes  a  phase  difference  between 
the  bending  and  torsional  oscillations  which  is  equal  to  90°  when  the  mass 
and  elastic  axes  are  coincident »  This  phase  difference  decreases  when  the 
mass  and  elastic  axes  are  not  coincident,  but  remains  substantial  in  the 
case  of  a  natural  mode  of  the  model  consisting  primarily  of  bending. 

(3)  Nonlinear  effects  for  large  motions  tend  to  change  the  natural  fre¬ 
quencies  of  the  system  slightly  and  introduce  some  coupling  between  the 
natural  vibration  modes  associated  with  solution  of  the  linearized  equations. 
A  limited  amount  of  results  did  not  provide  any  evidence  of  decay  or  diver¬ 
gence  of  the  free  vibrations  of  the  model . 
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APPENDIX  A 


DIFFERENTIAL  EQUATIONS  OF  MOTION 


The  differential  equations  for  free  motion  of  a  rotating  twisted  blade 
with  offset  mass  and  elastic  axes  are,  from  Ref.  2,  with  some  changes  in 
notation, 


-{[ 


} 


[GJe+TkA+EB^p'  f  ]</>'  -  EB2P' (6y  cos  P+5”  sin  0) 


+  TeA(5y  sin  P-5Z  cos  3)  +  ^pxe(-5y  sin  p+5^  cos  p) 

+  tf^peCsin  3)5y  +n  2p[(k|-k2)  cos  2(3+ee0  cos  p]^ 

2  2’)"  ”  2  1 
+  pCk^+k^p  -  pe(6y  sin  P-5Z  cos  P)  =  +  (TkAP' ) 

-  fl2p [(k|-k2)sin  P  cos  P+ee0  sin  p] 


((Eli  cos2  P+EI2  sin2  p)&£  +  (EI2-EIi)sin  p(cos  p)s£ 

-  TeAjzf  cos  p  -  EBaP'^'  sin  pj  "  -  (T&£)'  -  (ft2pxejz(  cos  p) ' 

+  p(&z+e^  cos  P)  =  (TeA  sin  p)"  +  (ft  pxe  sin  p)  (Al) 

••  •*  II  o  i 

+  p(5z+ep  cos  p)  =  (TeA  sin  p)  +  (ft  pxe  sin  p) 

QEI2-EI1)sin  p( cos  p)8£  +  (El!  sin2  P+EI2  cos2  p)5y 
+  TeA^  sin  p  -  EB2p'^'  cos  P)J  -  (T5y)’  +  (f^pxej^  sin  P) 

+  sfp erf  sin  P  +  p (Sy-eji?  sin  P)  -  tfpby  =  +(TeA  cos  p) 

+  (fl2pxe  cos  p)'  +  J^pCeQ+e  cos  P) 


An  explanation  of  the  origin  of  the  various  terms  in  the  equations  is 
given  in  Ref.  2.  The  integrals  which  define  the  section  constants  Bi  and  B2 


6h 


are  given  below 


Bi  =  ^'(T]2+C2-k|)(T12+^)dA 
B2  =  £  (n2+^2-k2)T)<lA. 


(A2) 


All  other  symbols  are  defined  in  the  list  of  symbols.  The  coordinate  system 
is  as  shown  in  Fig.  2.1. 

It  should  be  noted  that  Eq_s .  (Al)  are  for  small  displacements  from  the 
undeformed  configuration  of  the  blade  when  it  is  not  rotating.  The  analysis 
of  the  present  report  linearizes  the  problem  with  respect  to  small  displace¬ 
ments  from  the  steady-state  deformed  configuration  of  the  rotating  blade. 
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APPENDIX  B 


RESULTANT  LOADINGS 

The  resultant  loadings  per  unit  length  in  the  x,  £,  and  tj  directions  have 
been  obtained  in  Ref.  2  for  a  rotating  twisted  blade  with  offset  mass  and 
elastic  axes.  The  loads  include  the  inertial,  centrifugal  and  Coriolis  force 
terms.  In  the  notation  of  the  present  report  they  are 

px  =  -p(u-^u)  -  2HpSi  sin  p  +  2ftp$2  cos  p 

-  peSiP1  +  pe52  +  S^peBip*  -  J^peB^  -  20pe$  sin  p  +  ft2  px 

•• 

Ptj  =  -p52  +  ^pC-Bi  sin  p  cos  p+62  cos2  p+eQ  cos  p) 

-  2flpu  cos  p  +  ^pe  cos2  p  -  f^pe^  sin  p  cos  p 

-  fl^pe^  sin  P  -  f^pe^  sin  p  cos  p  +  2flpe  cos  p^BiP’+B^ 

p^  =  -pBi  -  I32p(-6i  sin2  p+S2  sin  pcosP+e0sin  p)  +  2flpu  sin  P 

•  • 

-  pejrf  -  ft2 pe  sin  p  cos  p  +  J^pe^  sin2  p  -  2fipe  sin  P(-5iP’+§2) 

(Bl) 

-  & peQ^  cos  P  -  fi2 pe^  cos2  p 

q.x  =  -A2 pe  [(-5i  sin  P+62  cos  p+eo)  sin  P  +  eQ^  cos  pj 

+  pe(-5i+20u  sin  P)  -  ST2  (jlg-I^)  sin  p  cos  p  +  cos  2p] 

“  sin  P-BiP1 cos  P+62  cos  p-62p'  sin  p)sin  P  cos  p 

+  2fi(l^  sin2  P+Ifj  cos2  P)(6i  cos  P-61P'  sin  p+82  sin  P+&2P'  cos  p) 

It,  =  n2IT)(5i+52B' )  -  IT)(Bi+62P')  +  2fll  $  COS  P 
q.£  =  -J^peCx+u)  +  peii  -  2pefl(-&1  sin  P+62  cos  P) 

-  flPl^BiP'  +  tiFlf&z  +  IjBiP’  t  I582  +  2«I^  sin  P 
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In  the  following,  we  eliminate  terms  in  px  which  are  dependent  on  dis¬ 
placement  variables  and  their  derivatives,  since  these  lead  to  nonlinearities 
in  subsequent  analysis,  and  terms  in  all  force  and  moment  expressions  which 
involve  u  and  its  derivatives  and  first  derivatives  of  the  remaining  displace¬ 
ment  variables.  We  also  eliminate  terms  involving  P',  since  these  terms  arise 
when  Si  and  &a  are  derivatives  referred  to  the  axes  q  and  t,  rotating  about  the 
x  axis.  In  the  lumped  parameter  treatment,  Si  and  62  can  be  considered  to  be 
derivatives  with  respect  to  locally  fixed  axes. 


Px 


=  sf  px 


p  =  -p&2  -  cFp  sin  p  (cos  p&i  +  ft  p  (cos2  pS2 

T) 

-  ft2p  sin  P(2e  cos  P+e0)^  +  ^P  cos  p(e  cos  P+e0) 

p^  =  -P&1  +  (sin2  "  fi2P  sin  P  fc°s  $  62  “  pe^ 

+  fi2p'jje(sin2  p-cos2  p)  -  eQ  cos  pj-  -  fi2p  sin  p(e  cos  P+e0) 

=  -pe&'i  +  Cfpe  (sin2  pBi  -  J^pe  sin  p  (cos  $02 

-  ft2  "jpee0  cos  p+(l£-I^)  cos  2p^J-jz(  -  (I£+It])^ 

-  ft2  sin  P'|'pee0+(l^-lT|)  cos  pj' 


qn  =  fl2!^  -  1^6 1 


q^  =  I  ^&2  “  I  ^&2  *■  ^  pex 


Now,  if  we  consider  the  matrix  equation 


(B2) 


=  (F  Jn^n 


(B3) 


where  (A)  is 


<<o 


(M  =< 


rv? 


Mi 

81 

61 
v2 

M2 

62 
52 
Q 

t 


V 


(B4) 


and.  (A}n  refers  to^the  value  of  these  quantities  at  station  N  just  outboard 
of  the  mass,  and  C ^ ) jj  refers  to  the  values  just  inboard  of  the  mass,  it  fol- 
lows  that 


Vi<N)  =  Vi(N)  +  p^N)i(N) 

=  VX(N)  +  p^M^ocP+fi2  sln2  p)51^N)  .  pWt(*)&  cos  p  £in  py8JN) 

+  pW/W  je^N^  (ac?+(sin2  f3-cos2  p)ft2}+e0ft2  cos  p  ]  *«(N) 

-  p(NMW)ft2  sin  3(eQ(N)+e(N)  cos  3) 

».<">  -  m,'”1  * »("VB)  -  P  <B  W»>  <B; 

TJ  A 

-  -  c(B><<BWBWB)  , 

W>  .  6;<"> 
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v«(H)  =  v2(n)  +  p 

=  +V2<N)  -p^  sinMeos 
+  p (NMN)(tu2+J22  cos2  p)S2^ 

-  p(H)|(H)q2  sin  p(e^2  cos 

*  cos  p{e0(^)+e<^>  cos  f» 


i'"' 

+ 

11 

(H)  (H)  (H)  (H)  (H) 

<1^  l  -  px  i  02 

.  Mp'"'  - 

00, (N)  (N)  a  (N)  (N)  (M)(H)  (HLa 

p  i  x  ii  62  “pzxeir^ 

.  6s<« 

(B5 ) 

4,»> 

CM 

cO 

II 

A(n) 

ft(H)  A 

(H).(H) 

Q 

=  Q  + 

1 

sin2  0)6i^ 

-  »W,W,I‘V  cos  e  jsin  etei"’  *  Q(,,) 

+  [(Vi5)<»)«ttV  ♦  (Vi;;)®*  (cos2  P-  sin2  p)l^fl? 

-  p(N)e(N)1(N)eo(N)fi2  cos  p]y(l) 

sin  p  cos  p  -  p^M^e^eo^n2  sin  p 

^M)  =  /*) 
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APPENDIX  C 


DEFORMATION  OF  A  BLADE  SEGMENT 


Consider  a  segment  of  a  weightless  beam  for  which  the  values  of  the  mo¬ 
ment  ,  torque,  shear,  and  tension  at  a  station  N  are  given  as  Mi^),  M2^> 

Q  N  >  Vi'  ,  V2'N'  and  T^.  Then  moments  and  torques  at  other  points  along 
the  segment  are 


Mi  =  Mi^)  +  Vjl^s  +  T^Bi  U(^V 

M2  =  m2W  +  v2^k)s  +  t(n)62  +  w(nV  (Cl) 


(n)  (n) 

where  Uv  '  and  Wv  '  are  the  contributions  of  the  centrifugal  force  coupling 
as  explained  in  Appendix  D,  s  is  the  longitudinal  coordinate  measured  from 
station  N  toward  the  root,  and  primes  denote  differentiation  with  respect 
to  s.  These  same  quantities  in  terms  of  5i  and  62  for  a  twisted  blade  are 


Mi  =  )2&i] 

M2  =  EI2^2-2p'&i-p"&i-0')282j  -  TeA  -  EB2pV  (C2) 

Q  =  -  [GJe-WkS+EBiO')2]  (*!'  -  Tk2p' 

+  EB2p’  [s2-2p,5i-p"8i-(p' )28l| 


For  a  straight  segment  Eqs>.  (C2)  reduce  to 


Mi  =  EI1S1 

M2  =  EI2&2  -  TeA  -  EB2P  V  (C3) 

Q  =  -  [GJe+Tk2+EB!(p’)2]  <!'  -  Tkjp’  +EB2(p’)62 
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Eliminating  steady-state  terms  and  combining  Eqs.(Cl)  and  (C3), 

“  a2^  =  bis  +  b2 

&2  “  a3^2  “  a4^ 1  -  a5$z(  =  bgS  +  b4 
a6&l  -  a7&2  +  a8^’  =  -t>5 


where 


ai  = 


a2 


a3 


a4  — 


t(k) 

Eli 

U« 

Eli 

t(n) 

ei2 

EBgp * 

EI2 


a  5  =  * 


W 


(N) 


a6 


EI2 

(N) 


U 


bi 


132 


133 


1>4 


Y.W 

El! 


Eli 


ei2 


Ela 


Ds 


ay  =  EB2  P 5 

a8  =  A(N)  =  GJe  +  Tk®  +  EBi(P')2 


(C4) 


These  equations  may  "be  solved  as  they  stand;  however,  considerable  sim¬ 
plification  can  he  achieved  at  a  modest  sacrifice  in  accuracy  by  replacing 
the  term  a5jzf  by  asjz^11).  This  modifies  the  second  of  Eqs.  (C4)  to  read 


82  -  a352  -  a4^'  =  b3s  +  be 


6  (C5) 
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where 


be  =  +^1 
ei2  ei2 

The  characteristic  equation  for  the  new  set  "becomes 

f0P4  +  ^iP2  +  *2  =0 


where 


fo  =  (a8-a4a7) 


fi  =  a2a6  -  ae(ai+a3)  +  aia4a7 


f2  =  0-10*38.8  -  o.2a3ae 


(06) 


The  roots  of  this  equation  are  then 


2 

Pi, 2 


(C7) 


(  JJ  ) 

When  the  solution  is  carried  out,  with  the  elements  of  (A)  '  as  initial  values, 

and  s  is  taken  equal  to  i,  the  following  form  for  the  solution  results 


bi(N+1) 


(N+i) 

Ol 


I (N+l ) 
&2 


E3jAj 


I 

1=1 


E7JAj 


(c8) 
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.  ,  10 
(N+l)  V1 
52  -  2^  EajAj 

j=l 

10 

^N+1>  =  £  Wj 
J-l 


(c8) 


where  the  Eij  are  presented  in  Eq.  (2.6)  in  the  body  of  the  report  and  Aj  are 
the  elements  of  (A}(N)  given  in  Eq.  (2.1). 

The  remaining  E^j  elements  are  found  from  the  following  equations  which 
apply  across  each  hay 


Vi0M) 

.  Vi®) 

.  Ml®>  *  V.®1/ 

v2®+1> 

II 

3 

(C9) 

=  M2®^  +  v»®>i 

,®w> 

-  ,®> 

' 

It  should  he  noted  that  the  bending  moment  and  torque  quantities  in  Eqs. 
(C9)  hy  definition  do  not  include  the  contributions  of  the  centrifugal  force 
displaced  to  the  elastic  axis  and  of  centrifugal  coupling.  These  contribu¬ 
tions  are  introduced  separately  through  satisfaction  of  Eqs.  (Cl)  in  the  solu¬ 
tion  for  the  deformation  variables. 
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APPENDIX  D 


CENTRIFUGAL  FORCE  COUPLING 


As  shown  in  Ref.  2 ,  there  is  a  type  of  coupling  between  bending  and  tor¬ 
sion  associated  with  the  presence  of  centrifugal  forces.  Explicit  considera¬ 
tion  must  be  given  to  the  derivation  of  the  terms  associated  with  this  coupling. 

If  xi  is  used  to  denote  the  station  where  centrifugal  force  is  acting  and 
x  the  station  where  bending  moment  is  measured,  the  components  of  bending  mo¬ 
ment  associated  with  offset  of  mass  center  from  the  elastic  axis  of  the  ro¬ 
tating  blade  may  be  written  as  follows: 


Mi  =  -  cos(p+$0  fl2pixiei  sin(Pi+^i 


)dxi 


4c 
R 


sin(p+^)^  f^piXiei  cos(Pi+^i  )dxi 


Ms  =  -  sin (&+$)  A2 PxXiex  sin  (PWi)dxi 

R 

(0+?O  f  r^pixxei  cos(£+^i)dxi 


(Dl) 


-  cos 


4c 


where  subscript  1  refers  to  values  at  xi- 

Assuming  f  to  be  a  small  angle  and  eliminating  higher  order  terms,  Eqs . 
(Dl)  become, 

»  rE 

Mi  =  -  cos  p  /  ^piXiei  sin  Pidxi  +  $  sin  p  /  ^PiXiei  sin  Pidxx 

R 


-  cos  ^PiXiei^i  cos  Pidxi 

R 

ld 


R  R 

+  sin  P  f  ft2 piXiei  cos  Pidxi  +  l  cos  P  f  fl^piXiei  cos  Pidxi 


sin  0  /  ^pixiei^i  sin  Pidxi 


(D2) 
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M2  =  -  sin  P  j  fi2piXiei  sin  Pidxi  -  $  cos  P  J  ^PiXi^i  sin  Pidxi 


-  sin  P  /  P^piXiei0i  cos  Pidxx 


R 


cos  p  J  fi^pixiex  cos  Pidxi  +  ^  sin  p J  ^2PiXiei  cos  Pidxi 
R 

+  cos  p  r  sin  pidxi  (^) 


In  each  case,  the  first  and  fourth  terms  represent  steady-state  moments. 
An  examination  of  Appendix  B  shows  that  the  term  -f^pex  in  will  give  rise 
to  these  moments.  They  are  taken  into  account  through  the  element  d6  in  the 
{d}  matrix. 

The  third  and  sixth  terms  in  each  component  represent  the  effect. of  tor¬ 
sional  displacement  of  the  Blade  mass  on  "bending  about  the  torsionally  undis¬ 
placed  positions  of  the  T)  and  I  axes  in  the  Mi  and  M2  components  respectively. 
This  effect  can  he  taken  into  account  by  an  appropriate  modification  of  the 
[F]  matrix,  incorporating  a  change  in  bending  moment  across  each  mass  given 

*>y 


AMi  =  -  f^pixe^ 


(D3) 


yielding  the  element 


F210  —  -  £2 2 pixe 


(D4) 


The  second  and  fifth  terms  represent  the  effect  of  centrifugal  forces 
acting  on  the  torsionally  undisplaced  masses  between  x  and  R  on  bending  about 
the  torsionally  displaced  T)  and  £  axes  in  the  Mi  and  M2  components  respectively. 
They  must  be  taken  into  account  in  the  development  of  the  [E]  matrix.  In  terms 
of  the  lumped  mass  model,  the  contribution  to  the  bending  moments  in  the  bay 
between  the  nth  and  (n+l)th  masses  is 
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Mi  = 


^in  ^ 
'  i=: 


u<“V 


M2  =  [  -COS 


i± PiXiei  Sin  P-^+cos  Pn  )  ^  Pi*ixiei  cos  Pi  )  i 

1  i=l 


n  n  \ 

Pn  ^  ^PiijXjei  sin  Pi+sin  Pn  ^  ^PiiiXiei  cos  PjJ  ^ 
i=l  i=l  ' 


(D5) 


w(% 


There  is  correspondingly  an  effect  of  bending  on  torque.  It  is  associated 
with  the  fact  that  with  a  bending  slope  6^  at  station  x  the  centrifugal  force 
ftspixidxi  on  an  element  dxi  at  station  xi  outboard  of  station  x  has  a  component 
-f^pixidxi&i  normal  to  the  rj-axis  in  the  plane  of  the  cross  section  at  station 
x.  The  moment  arm  of  this  force  about  the  elastic  axis  is  ei  cos(p-Pi),  so 
that  the  contribution  to  the  torque  from  this  source  may  be  written, 


K 

Q  =  -  ^PiXiex  cos  (p-pi  )dxi 


R 


R 


-  [  sin  P 


J  n2PiXiei  sin  Pidxi+cos  p  /  P^PiXiex  cos  Pidxi)&i  (D6) 


9 


Applying  this  result  to  the  lumped  mass  model,  the  contribution  to  the 
torque  in  the  bay  between  the  nth  and  (n+l)th  masses  is 


/  n  n  \ 

Q  =  -  (sin  pn  ^  f^PiljX-je  i  sin  p  +cos  Pn^  ^2PiJ^ixiei  cos  P  y  Si 

\  i=l  i=l  ' 


(D7) 


and  is  taken  into  account  in  the  development  of  the  [E]  matrix  in  Appendix  C. 
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